AFRL-SN-WP-TR-1999-1131 


INVESTIGATION  OF  POLYMER 
OPTICAL  WAVEGUIDE  DEVICES 


JAMES  G.  GROTE 
JOSEPH  E.  BRANDELIK 


AFRL/SNDI 

SENSORS  DIRECTORATE 
AIR  FORCE  RESEARCH  LABORATORY 
2241  AVIONICS  CIRCLE,  STE.  20 
WRIGHT-PATTERSON  AFB,  OH  45433-7322 


NOVEMBER  1999 


FINAL  REPORT  FOR  05/05/1995  -  01/10/1999 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED 


SENSORS  DIRECTORATE 

AIR  FORCE  RESEARCH  LABORATORY 

AIR  FORCE  MATERIEL  COMMAND 


20000525  023 


WRIGHT-PATTERSON  AIR  FORCE  BASE  OH  45433-7318 


dtic  quality  nmmamB  i 


NOTICE 


USING  GOVERNMENT  DRAWINGS,  SPECIFICATIONS,  OR  OTHER  DATA 
INCLUDED  IN  THIS  DOCUMENT  FOR  ANY  PURPOSE  OTHER  THAN 
GOVERNMENT  PROCUREMENT  DOES  NOT  IN  ANY  WAY  OBLIGATE  THE  US 
GOVERNMENT.  THE  FACT  THAT  THE  GOVERNMENT  FORMULATED  OR 
SUPPLIED  THE  DRAWINGS,  SPECIFICATIONS,  OR  OTHER  DATA  DOES  NOT 
LICENSE  THE  HOLDER  OR  ANY  OTHER  PERSON  OR  CORPORATION;  OR 
CONVEY  ANY  RIGHTS  OR  PERMISSION  TO  MANUFACTURE,  USE,  OR  SELL 
ANY  PATENTED  INVENTION  THAT  MAY  RELATE  TO  THEM. 

THIS  REPORT  IS  RELEASABLE  TO  THE  NATIONAL  TECHNICAL 
INFORMATION  SERVICE  (NTIS).  AT  NTIS,  IT  WELL  BE  AVAILABLE  TO  THE 
GENERAL  PUBLIC,  INCLUDING  FOREIGN  NATIONS. 

THIS  TECHNICAL  REPORT  HAS  BEEN  REVIEWED  AND  IS  APPROVED 
FOR  PUBLICATION. 


JOSEPH  E.  BRANDELIK 
Multi-Chip  Integration  Branch 
Aerospace  Components  Division 


CLAN  J.  TEWKSBURY,  Chiei 
Multi-Chip  Integration  BrancI 
Aerospace  Components  Division 


ROBERT  T.  KEMERLEY,  Chief 
Aerospace  Components  Division 
Sensors  Directorate 


Do  not  return  copies  of  this  report  unless  contractual  obligations  or  notice  on  a  specific 
document  requires  its  return. 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing 
the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information 
Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  222024  302,  and  to  the  Office  of  Management  and  Budget,  Paperwork  Reduction  Project  (0704-0188),  Washington,  DC  20503. 


1.  AGENCY  USE  ONLY  (Leave blank)  2.  REPORT  DATE  3.  REPORT  TYPE  AND  DATES  COVERED 

NOVEMBER  1999  FINAL  REPORT  FOR  05/05/1995  -  01/10/1999 


3.  REPORT  TYPE  AND  DATES  COVERED 


4.  TITLE  AND  SUBTITLE  5.  FUNDING  NUMBERS 

INVESTIGATION  OF  POLYMER  OPTICAL  WAVEGUIDE  DEVICES  C  IN-HOUSE 

PE  62204 
PR  ARPE 

6.  AUTHOR(S)  TA  01 

JAMES  G.  GROTE  qj 

JOSEPH  E.BRANDELIK 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

AFRL/SNDI 

SENSORS  DIRECTORATE 

AIR  FORCE  RESEARCH  LABORATORY 

2241  AVIONICS  CIRCLE,  STE.  20 

WRIGHT-PATTERSON  AFB.  OH  45433-7322 _ 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

SENSORS  DIRECTORATE 
AIR  FORCE  RESEARCH  LABORATORY 
AIR  FORCE  MATERIEL  COMMAND 
WRIGHT-PATTERSON  AFB,  OH  45433-7318 

POC:  JOSEPH  E.  BRANDELIK.  AFRL/SNDI.  937-255-4557  EXT.  3446 


11.  SUPPLEMENTARY  NOTES 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


10.  SPONSORING/MONITORING 
AGENCY  REPORT  NUMBER 


AFRL-SN- WP-TR- 1 999- 1 131 


THIS  IN-HOUSE  EFFORT  WAS  FUNDED  IN  PART  BY  DARPA. 


12a.  DISTRIBUTION  AVAILABILITY  STATEMENT 


APPROVED  FOR  PUBLIC  RELEASE,  DISTRIBUTION  UNLIMITED. 


13.  ABSTRACT  (Maximum  200  words) 

This  report  is  an  initial  investigation  into  polymer  waveguides.  The  CAD  program  BEAMPROP  is  utilized  for  modeling  of 
waveguides.  Passive  multi-mode  and  single-mode  waveguides  are  considered.  For  modeling  nonlinear  optical  polymers,  in 
addition  to  the  BEAMPROP AGATION  (BEAMPROP)  method,  the  effective  index  method  and  the  MARCATILI  model  are 
developed.  It  is  shown  that  a  conducting  cladding  layer  for  the  rectangular  waveguides  could  significantly  improve  the 
operation  performance  of  nonlinear  optical,  poled  polymer  waveguides. 


14.  SUBJECT  TERMS 


Optical  Polymers,  Waveguides,  Waveguide  Modeling 


15.  NUMBER  OF  PAGES 


|16.  PRICE  CODE 


17.  SECURITY  CLASSIFICATION 
OF  REPORT 

UNCLASSIFIED 


18.  SECURITY  CLASSIFICATION 
OF  THIS  PAGE 

UNCLASSIFIED 


19.  SECURITY  CLASSIFICATION 
OF  ABSTRACT 

UNCLASSIFIED 


20.  LIMITATION  OF 
ABSTRACT 


Standard  Form  298  (Rev.  2-89)  (EG) 

Prescribed  by  ANSI  Std.  23S.18 

Designed  using  Perform  Pro,  WHS/DIOR,  Oct  94 


TABLE  OF  CONTENTS 


LIST  OF  FIGURES  v 

LIST  OF  TABLES  viii 

FOREWARD  ix 

ACKNOWLEDGMENTS  x 

SUMMARY  1 

1.  INTRODUCTION  2 

2.  THE  POINT  PROGRAM  5 

2. 1  Related  Program  6 

2.2  Program  Approach  7 

3.  THE  COMPUTER  PROGRAM  BEAMPROP  9 

3.1  Beamprop  Enhancements  10 

3.2  Theoretical  Background  1 1 

4.  BEAMPROP  RESULTS  14 

4.1  Calculation  Of  Crosstalk  In  The  Optical  Backplane  14 

4.2  Simulations  On  Multimode  S-Bends  16 

4.3  Simulations  For  Single-Mode  Directional  Couplers  18 

5.  DESIGN  AND  FABRICATION  OF  NON-LINEAR  OPTICAL  POLYMER 

INTEGRATED  OPTICAL  DEVICES  20 

5.1  Electro-Optic  Properties  Of  Nonlinear  Optical  Polymers  20 

5.1.1  Linear  Electro-Optic  Effect  22 

5.1.2  Index  Ellipsoid  26 

5.1.3  Electro-Optic  Coefficients  29 

5.1.4  Induced  Refractive  Indices  and  Birefringence  34 

5.2  Electro-Optic  Modulators  37 

5.3  Directional  Couplers  40 

5.3.1  Zero-Gap  Directional  Couplers  43 

5.3.2  Wavefront  Tilt  46 

5.3.3  Interaction  Length  Of  The  Directional  Coupler  47 

5.3.4  Rectangular  Dielectric  Waveguide  Model  48 

iii 


TABLE  OF  CONTENTS  (continued) 

5.3.5  Marcatili  TEpq  Modes  50 

5.3.6  Marcatili  TMpq  Modes  55 

5.3.7  Effective  Index  Method  56 

5.3.8  Effective  Index  Method  TEpq  Modes  58 

5.3.9  Effective  Index  Method  TMpq  Modes  64 

5.3.10  Beam  Propagation  Method  66 

5.3. 1 1  Refractive  Index  Of  The  Waveguide  Cladding  72 

5.4  Device  Fabrication  74 

5.4.1  Photobleach  74 

5.4.2  Etch  And  Fill  76 

5.4.3  Corona  Poling  76 

5.5  Practical  Nonlinear  Optical  Polymer  Electro-Optic  Devices  78 

5.5.1  Electrode  Separation  79 

5.5.2  Passive  Zero-Gap  Directional  Coupler  81 

5.5.3  Active  Zero-Gap  Directional  Coupler  84 

5.5.4  Nonlinear  Polymer  To  Passive  Polymer  Coupling  87 

6.  CONCLUSIONS  88 

7.  REFERENCES  89 


IV 


LIST  OF  FIGURES 


Figure  1 . 

Figure  2. 
Figure  3. 
Figure  4. 
Figure  5. 
Figure  6. 
Figure  7. 
Figure  8. 

Figure  9. 

Figure  10. 

Figure  1 1. 
Figure  12. 
Figure  13. 


Power  Requirements  for  Optical  and  Electrical  Interconnects  versus 
Data  Rate  [32]. 

Simulations  for  an  array  of  five  waveguides. 

Results  for  offsets  of  1, 2,  and  5  pm. 

Simulation  results  for  2.5'S-bends. 

Simulation  results  for  7.5°  S-bends. 

Simulation  results  for  smooth  1.5°  S-bends. 

Simulation  of  a  single-mode  directional  coupler. 

Index  Surfaces  for  (a)  a  Positive  Uniaxial  Crystal  and  (b)  a  Negative 
Uniaxial  Crystal  [29]. 

Inner  Ellipse  is  the  Intersection  of  Index  Ellipsoid  with  the  Plane 
Perpendicular  to  s  [15]. 

Transverse  Electro-Optic  Modulator  using  a  Tetragonal  4  Group 
Crystal. 

Operation  of  a  Mach  Zehnder  Interferometer. 

Operation  of  a  Directional  Coupler. 

Operation  of  a  Zero-Gap  Directional  Coupler. 


LIST  OF  FIGURES  (continued) 


Figure  14.  Rectangular  Waveguide  Directional  Coupler.  47 

Figure  15.  Rectangular  Waveguide  Directional  Coupler  Cross  Section.  49 

Figure  1 6.  Field  Intensity  Profiles  of  the  Symmetric  and  Antisymmetric  Modes 

in  a  Rectangular  Directional  Coupler.  52 

Figure  17.  Trajectory  of  a  Guided  Wave  in  a  Slab  Waveguide.  56 

Figure  1 8.  Analytical  Model  for  the  Effective  Index  Method  for  TEpq  Modes.  57 

Figure  19.  Amplitudes  of  the  Propagating  Fields  in  an  Integrated  Optic  S-Bend 

Waveguide.  72 

Figure  20.  Schematic  of  Triple  Stack  Including  Conductive  Polymer  Charge 

Sheets.  81 

Figure  2 1 .  Schematic  of  a  Passive  Zero-Gap  Directional  Coupler  With  an  NLO 

Polymer  Interaction  Region  and  Passive  Polymer  Input/Output 
Channels.  82 

Figure  22.  Beam  Propagation  Model  of  Passive  Zero-Gap  Directional  Coupler 

Using  NLO  and  Passive  Polymers  With  Matching  Indices.  83 

Figure  23.  Beam  Propagation  Model  of  Passive  Zero-Gap  Directional  Coupler 

Using  NLO  and  Passive  Polymers  With  0.01  Mismatched  Indices.  84 


VI 


LIST  OF  FIGURES  (continued) 


Figure  24. 


Figure  25. 


Beam  Propagation  Model  Of  Active  Zero-Gap  Directional  Coupler 
Using  NLO  And  Passive  Polymers  With  Matching  Indices.  85 


Beam  Propagation  Model  Of  Active  Zero-Gap  Directional  Coupler 
Using  NLO  And  Passive  Polymers  With  0.01  Mismatched  Indices.  86 


vii 


LIST  OF  TABLES 


Table  1. 

Major  Programs  On  High  Speed  Parallel  Optical  Data  Links. 

8 

Table  2. 

NLO  Materials  Comparison. 

79 

Table  3. 

NLO  Polymer  Comparison. 

80 

viii 


FOREWORD 


While  passive  polymer  waveguides  are  addressed  in  this  manuscript,  this  report 
focuses  on  the  modeling  and  fabrication  of  practical  and  competitive  non-linear  optic 
(NLO)  polymer  integrated  optic  devices  for  computer  based  applications.  The  devices 
covered  will  include  mach  Zehnder  interferometers,  conventional  directional  couplers 
and  zero-gap  directional  couplers.  Modeling  techniques  will  include  dielectric 
waveguide,  effective  index  and  beam  propagation  methods.  Fabrication  techniques, 
including  photobleach,  and  etch  and  fill,  will  be  detailed.  State-of-the-art  design  schemes 
exhibiting  the  potential  to  lower  the  switching  voltage  of  NLO  polymer  integrated  optic 
devices  will  be  investigated,  including  the  use  of  conductive  polymers.  The  NLO 
polymers  presented  in  this  manuscript  will  be  those  readily  available.  Coupling 
techniques  used  to  achieve  good  coupling  efficiency  between  passive  polymers,  utilized 
for  general  routing,  and  NLO  polymers,  used  for  active  switching,  will  also  be  addressed. 


IX 


ACKNOWLEDGMENTS 


The  optically  passive  work  is  based  in  part  on  discussions  with  Yung  Liu  of  General 
Electric  Corporate  Research  and  Development  (GECRD),  Mario  Ghezzo  of  GECRD 
Matthew  Nielsen  of  GECRD  and  Antonio  Crespo  of  the  Air  Force  Research 
Laboratory/Sensors  Directorate  (AFRL/SNDI).  The  optically  active  polymer  work  is  based 
in  part  on  discussions  with  Julian  Bristow  of  Honeywell,  Stephen  Caracci  of  Air  Force 
Research  Laboratory/Materials  and  Manufacturing  Directorate  (AFRL/MLPO),  Bruce 
Reinhardt  of  AFRL/MLBP,  Mohammad  Karim  of  the  University  of  Dayton  (UD),  Vincent 
Dominic  of  UD,  Charles  Lee  of  the  Air  Force  Research  Laboratory/Air  Force  Office  of 
Scientific  Research  (AFRL/AFOSR),  Ray  Chen  of  the  University  of  Texas  at  Austin,  Larry 
Dalton  of  the  University  of  Southern  California,  Paul  Ashley  of  the  US  Army  Missile 
Command,  and  Charles  Sullivan  of  Sandia  National  Laboratories  Their  input  was  both 
welcomed  and  invaluable.  Sincere  appreciation  is  extended  to  the  members  of  AFRL/SN 
for  their  help  and  support  in  device  fabrication  and  test  equipment  design  and  fabrication. 
Members  of  Technology  Scientific  Services  and  Wright  State  University,  working  for 
AFRL/SN  also  provided  assistance  in  device  fabrication,  generation  of  waveguide  masks 
and  design  and  fabrication  of  the  bases  and  mounts  for  the  test  equipment. 

This  manuscript  covers  in-house  work  spanning  a  number  of  years  and  was 
funded  by  the  Air  Force  Research  Laboratory  (AFRL)  and  the  Defense  Advanced 
Research  Projects  Agency  (DARPA). 


SUMMARY 


This  document  covers  optically  passive  and  optically  active  polymer  waveguides. 
The  passive  polymers  that  are  modeled  are  multimode  and  the  work  is  an  outgrowth  of 
the  General  Electric  POINT  program  (Air  Force  contract  F33615-94-C-1530).  Modeling 
of  the  waveguides  was  through  the  use  of  the  BeamPROP  CAD  program.  Information 
transfer  using  passive  polymers  requires  modulation  of  a  laser.  Optically  active 
waveguides  can  provide  their  own  modulation.  The  modeling  and  fabrication  of  non¬ 
linear  optic  (NLO)  polymer  integrated  optic  devices  for  computer  based  applications  are 
presented.  The  devices  covered  will  include  Mach-Zehnder  interferometers, 
conventional  directional  couplers  and  zero-gap  directional  couplers.  Modeling 
techniques  include  the  dielectric  waveguide,  effective  index  and  beam  propagation 
methods.  Fabrication  techniques,  including  photobleach,  and  etch  and  fill,  are  detailed. 
Design  schemes  that  exhibit  the  potential  to  lower  the  switching  voltage  of  NLO  polymer 
integrated  optic  devices  are  investigated,  including  the  use  of  conductive  polymers. 
Coupling  techniques  used  to  achieve  good  coupling  efficiency  between  passive  polymers, 
utilized  for  general  routing,  and  NLO  polymers,  used  for  active  switching,  are  also  be 
addressed. 
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1.0  INTRODUCTION 


As  electronic  and  optical  data  processing  circuits  approach  multi-gigahertz  operation 
rates,  the  need  for  multi-gigahertz  communication  between  these  processing  circuits 
becomes  essential.  This  communication  includes  chip-to-chip  on  a  common  multichip 
module  (MCM),  MCM  to  MCM  on  a  common  board,  and  board-to-board  through  a 
common  backplane.  The  data  rates  at  which  communication  through  electrical 
interconnects  can  be  transmitted  is  limited  to  the  low  gigahertz  rates,  unless  one  is  willing  to 
expend  large  amounts  of  power.  Electrical  wires  and  striplines  posses  frequency  dependent 
dissipative  losses  due  to  conductive  and  dielectric  properties,  and  frequency  dependent 
parasitic  coupling  losses  due  to  radiative  and  surface  wave  propagation  properties  [1].  The 
higher  the  frequency,  the  higher  the  losses.  This  is  where  optics  can  help  provide  solutions. 
Optical  interconnects  can  transmit  data  at  optical  frequencies  with  orders  of  magnitude 
lower  losses  than  their  electrical  counterparts.  The  power  requirements  for  optical  and 
electrical  interconnects  versus  data  rate  are  illustrated  in  Figure  1 . 


[32], 


The  criterion  used  to  generate  Fig.  1  were  as  follows:  waveguide  lengths  for  both 
the  electronic  and  optical  interconnects  were  50  mm.  The  total  number  of  interconnects  was 
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600  for  both  types.  For  electrical  interconnects  on  an  alumina  substrate,  waveguide  widths 
of  25  jim  and  waveguide  heights  of  650  |im  were  assumed.  The  effective  and  relative 
permittivities  of  the  microstrip  used  were  6.7  and  9.9,  respectively  [1].  The  capacitance 
used  per  line  was  6  pF  and  an  impedance  of  50  Q.  was  assumed.  The  power  versus  data  rate 
for  the  electrical  interconnects  was  determined  with  methods  given  by  Edwards  [1].  For 
optical  interconnects  on  a  semiconductor  substrate,  waveguide  widths  of  3  pm  and 
waveguide  heights  of  1  pm  were  assumed.  Waveguide  propagation  loss  of  0.2  dB/cm  was 
assumed  for  1  GHz  growing  as  0.08  dB/GHz  [2, 3, 4, 5, 6, 7, 8].  The  modulator  loss  was 
5  x  10'5  W/channel  for  1  GHz  growing  as  5  x  10‘5  W/GHz  [3].  The  detector  capacitance 
was  0.3  pF  and  the  detector  amplifier  capacitance  used  was  1.7  pF  [3].  The  power,  per 
detector,  was  1  x  10'3  W/channel  for  1  GHz  growing  as  1  x  10‘3  W/GHz  [3].  The  laser 
power  was  100  mW  for  1  GHz  decreasing  at  the  same  rate  as  the  waveguide  loss  was 
growing  per  GHz  [3].  Finally,  20%  electrical  efficiency  and  50%  optical  coupling 
efficiency  were  assumed  for  the  laser.  The  power  versus  data  rate  was  then  determined  by 
adding  up  the  power  usage  of  the  components. 

In  addition  to  the  lower  power  consumption  and  higher  bandwidth,  optical 
interconnects  offer  isolation  from  electro-magnetic  interference  (EMI)  [9].  Optical 
interconnects  are  starting  to  be  utilized  in  chassis-to-chassis  and  board-to-board 
communication  using  glass  optical  fibers  [2, 5, 9, 10, 11, 12].  Polymeric  optical  waveguides 
are  in  the  early  stages  of  application  to  interconnects.  On  the  smaller  scale  of  multichip 
modules  and  within-module,  integrated  optics  is  just  beginning  to  emerge  [2,  3, 5, 6,  8-40]. 
Optics  offers  the  additional  advantages  of  smaller  interconnect  widths  and  in-plane 
crossover  capability  with  very  low  crosstalk.  This  eliminates  the  need  for  air  bridges.  Less 
wafer  and  board  real  estate  will  be  required  for  optical  interconnects,  which  leaves  room  for 
additional  electronic  circuits  and  devices.  Fabrication  processes  for  optical  interconnects 
must  be  compatible  with  those  processes  used  for  manufacturing  electronic  interconnects 
and  devices.  This  is  important  from  the  standpoint  of  practicality  and  of  acceptability. 
Polymers  are  attractive  materials  because  they  are  currently  being  used  for  the  buffer  layers 
between  signal  planes  in  IC's  and  MCM's.  Nonlinear  optical  (NLO)  polymers  possess 
electro-optical  effects  which  can  be  utilized  to  fabricate  electro-optic  (EO)  switches. 
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Optical  switches  can  be  incorporated  with  optical  waveguides  and  other  EO  switches 
to  realize  digital  logic  operations  reconfigurable  interconnects.  Current  research  in 
developing  processing  techniques  for  polymer  optical  waveguides  and  devices  has  shown 
that  the  processing  techniques  for  these  optical  circuits  are  compatible  with  those  techniques 
used  for  the  fabrication  of  electronic  circuitry  [2, 3, 4, 5, 6, 7, 8, 21, 22, 35, 36, 39-65],  By 
using  NLO  polymers,  one  can  conceivably  house  the  entire  electronic/electro-optic  circuit 
on  a  single  substrate.  With  the  incorporation  of  NLO  polymer  electro-optic  circuits  and 
polymer  optical  waveguides,  the  multichip  module  or  board,  currently  entirely  electronic, 
could  perform  more  efficiently  with  less  power  consumption,  move  data  around  faster  with 
less  noise,  and  use  less  real  estate. 

There  is  also  the  potential  to  reconfigure  the  multichip  modules  to  perform  entirely 
different  operations,  or  to  route  signals  around  inoperative  or  defective  portions  of  the 
module  [2, 3,  8-1 1, 14, 20-22,  32, 35-39, 41, 52,  53]. 

Nonlinear  optic  polymers  have  several  attractive  characteristics  which  many 
researchers  have  tried  to  capitalize  on  over  the  past  decade.  These  include  a  high  electro¬ 
optic  coefficient  for  low  voltage  operation,  a  low  dielectric  constant  for  high  speed 
modulation,  low  temperature  processing  for  integration  of  optics  with  electronics,  excellent 
index  match  to  optical  fibers  and  the  potential  to  simplify  fabrication  for  lower  cost  [3, 4, 7, 
8, 21 , 22, 35, 36,  38-40, 42, 45-47, 49-70].  Current  properties  of  readily  available  NLO 
polymers  include  an  electro-optic  coefficient  r33  =  26  pm/V  at  X  =  830  nm,  a  propagation 
loss  of  1.5  dB/cm  at  X  =  830  nm  and  a  long  term  poling  retention  temperature  stability  of 
80°C  [3, 4, 7,  8,  21, 22, 35,  36,  38-40, 42, 45-47, 49-70].  These  properties  should  be 
adequate  for  demonstration  of  a  competitive  electro-optic  device. 
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2.0  The  Point  Program 


There  have  been  significant  advances  in  optical  materials,  devices,  and  packaging 
technologies,  making  optical  interconnects  the  viable  wide  band  interconnect  technology  of 
choice  [71].  These  enabling  technologies  are: 

1)  Vertical  cavity  surface  emitting  lasers  (VCSEL)  devices  and  low  cost  VCSEL  packaging 

2)  Polymer  based  optical  waveguide  and  fiber  materials 

3)  Low  cost  micro-optical  components,  connectors,  and  fabrication  processes. 

Some  of  these  new  optoelectronic  materials  and  devices  can  be  handled,  processed,  and 
packaged  much  like  IC  devices  using  planar  fabrication  and  batch  processing,  meeting  the  key 
requirements  for  low  cost,  large  volume  manufacturing. 

The  POINT  program  was  organized  in  response  to  DARPA’s  BAA-94-46,  "Affordable 
Optoelectronic  Module  Technology,”  and  was  under  contract  on  December  1,  1994.  The 
Program  kickoff  meeting  was  held  on  December  15,  1994  at  DARPA,  and  was  attended  by  all 
Principal  Investigators  -  General  Electric,  Honeywell,  AlliedSignal,  AMP,  Columbia  University, 
and  the  University  of  California  at  San  Diego  (UCSD).  The  goal  of  the  POINT  program  was  to 
develop  affordable  optoelectronic  packaging  technologies  utilizing  optical  polymer  media  for 
board  and  backplane  optical  interconnect  applications  in  military  and  commercial  systems.  The 
program's  technical  approach  was  to  leverage  electronic  packaging  and  fabrication  processes  - 
such  as  batch  process,  planar  fabrication  and  passive  alignment  -  to  significantly  reduce  the  cost 
of  optoelectronic  packaging.  The  program  development  focused  on: 

1. )  VCSEL  devices  and  the  related  packaging  technology  that  leverages  the  electronic  planar 

fabrication  process 

2. )  Passive  alignment  techniques  that  could  reduce  significantly  the  recurrence  cost  in  O/E 

packaging 
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3. )  Planar  polymer  waveguides  that  could  be  used  for  board  and  backplane  applications  for 

high  density 

4. )  High  speed  optical  interconnects  to  remove  the  interconnect  bottlenecks  among  multiple 

processors  in  a  computer  system. 

2.1  Related  Programs 

The  Optoelectronic  Technology  Consortium  (OETC)  [72]  was  a  DARPA  sponsored 
industrial  consortium.  Originally  formed  in  1992  by  GE  (later  transferred  to  Lockheed  Martin), 
ATT,  IBM,  and  Honeywell,  OETC  aimed  at  developing  key  optical  interconnect  technology  for 
computer  applications.  This  represented  the  first  major  R&D  effort  to  develop  32-channel 
parallel  optical  links  operating  at  500  Mbps  per  channel,  with  a  BER  <  1  X  10-14  and  average 
power  of  0.25  W/channel  dissipated  power.  The  link  included  32-channel  AlGaAs  VCSEL 
arrays,  driver  arrays,  GaAs  MESFET  receiver  arrays,  and  fiber  ribbon  arrays  with  MAC  II 
connectors.  This  program  is  credited  with  recognizing  the  significance  of  VCSEL  devices,  and  it 
accelerated  the  development  of  VCSEL  technology  from  a  laboratory  novelty  at  the  beginning  of 
the  program,  to  a  critical  optical  component  by  the  time  the  program  ended  in  1995,  meeting  all 
requirements  for  short-haul  optical  communication. 

The  POLO  program  [73]  consisted  of  BP,  AMP,  DuPont,  Spectra-Diode  Labs  (SDL),  and 
University  of  Southern  California  (USC),  and  was  sponsored  jointly  by  industry  and  DARPA. 
Started  in  1994,  it  aimed  to  reduce  the  cost  of  optoelectronic  module  technology.  POLO  led  to 
the  development  of  bi-directional  10-channel  parallel  optical  links  operating  at  1  Gbps  per 
channel. 

The  JITNEY  program  [74],  a  joint  development  effort  by  IBM  and  3M,  was  a  NIST- 
sponsored  ATP  program  to  develop  low  cost,  parallel  optical  modules  operating  at  20  channels 
and  1  Gbps/  channel.  JITNEY  used  20-channel  VCSEL  arrays,  CMOS  driver  ICs,  and  a  20- 
channel  GaAs  MESFET  receiver  array  OEIC.  The  modules  were  made  of  plastic  molded  optical 
couplers  to  interface  between  VCSEL  and  large  core  fibers. 
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OPTOBUS  is  a  commercial  product  developed  by  Motorola.  It  is  available  as  bi¬ 
directional  10-channel  transceiver  modules  operated  at  200  Mbps/channel. 

2.2  Program  Approach 

Optical  interconnect  technologies  can  be  classified  broadly  into  three  categories 
according  to  the  transmission,  data  rate,  and  distance: 

1 . )  For  long-haul  communication  over  several  kilometers  (such  as  telecommunication, 

wide  area  network,  single-mode  optic  fibers  and  long  wavelength),  using  InGaAsP 
edge-emitting  lasers  at  1.3  at  1.5  Jim.  Single-channel  point-to-point  interconnects  are 
typical. 

2. )  For  a  communication  distance  ranging  between  tens  and  hundreds  of  meters  (such  as 

local  area  network  interconnects  between  high  speed  servers,  workstation  clusters,  and 
supercomputers,  using  multimode  optic  fibers  and  short  wavelength).  AlGaAs/GaAs 
semiconductor  lasers  at  0.85  Jim  are  commonly  employed,  using  either  single  or 
multiple  channels.  The  OETC,  POLO,  OP'TOBUS,  and  JITNEY  programs  developed 
the  technologies  for  these  applications. 

3. )  For  short  distance  communication  less  than  a  few  meters  (such  as  within  a  cabinet,  a 

backplane,  or  an  interconnect  between  frames  and  boards),  employing  multichannel 
VCSEL  arrays  at  0.85  Jim  with  multimode  fibers  or  polymer  waveguides. 

The  POINT  program  approach  was  to  apply  emerging  enabling  technologies  with  the  intent 
of  developing  an  affordable  optical  interconnect  and  packaging  technology.  In  the  POINT 
program,  Honeywell's  state-of-the-art  VCSEL  devices  were  used  as  the  laser  source  for  the 
development  of  the  parallel  optical  interconnect  links.  The  VCSEL  is  the  single  most  important 
active  optical  component  developed  for  interconnect  applications  in  short-haul  data 
communication.  The  POINT  was  aimed  at  the  development  of  an  affordable  optoelectronic 
packaging  and  interconnect  technologies  using  optical  polymers  for  board  and  backplane 
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applications.  See  Table  1  for  a  comparison  of  programs.  The  expression,  BER  = 
0.5*erfc(SNR/(2*sqrt(2))),  was  used  to  determine  the  bit  error  rate  from  the  signal  to  noise  ratio. 


Table  1.  Major  Programs  On  High  Speed  Parallel  Optical  Data  Links. 


Project 

OETC 

OptoBus 

POLO 

Jitney 

POINT 

Sponsor 

DARPA 

Motorola 

DARPA/TRP 

NIST/ATP 

DARPA 

Team 

GE  (LMC), 

ATT,  IBM, 
Honeywell 

Motorola 

HP,  AMP,  SDL, 

use 

IBM,  3M 

GE,  Honeywell, 
Allied  Signal, 
AMP,  Columbia, 
UCSD 

#  channels 

32 

10/10 

10/10 

20 

12/12 

Data  rate  per 
channel  (Mbps) 

>500 

150 

625,1000 

500 

>350 

Total  data  rate 
(Gbps) 

10020 

1.5 

6 

10 

4.2 

Max  link  length 
(m) 

100 

30 

1-300 

<30 

Board, 

backplanecO.l 

Size  (mm) 

33x44x5 

38x38x13 

40x40 

28x18 

N?A 

Power  (W) 

Tx+Rx 

10 

1.7 

<3 

<3 

2 

B?  1 

140 

250 

250 

500 

250/100 

Lead  frame 

Lead  frame 

Lead  frame  SM 

User  defined 

MACH 

MT 

MT 

- 

MT 

VCSEL  (850) 

VCSEL  (850) 

VCSEL  (850) 

VCSEL  (850) 

VCSEL  (850) 

Fiber  coupling 

Fiber/SiOB  45° 
mirror 

1 

Moulded  Plastic 
couplers 

Polyguide  45° 
mirror 

Butt  coupled 

Polyguide  & 
polymer 
waveguide  45° 
mirror 

Fiber/pitch 
dimensions  (pm) 

62.5/140 

62.5 

62.5/125 

200 

50/100 

(rectangular 

waveguides) 

Electrical 

interface 

LVDS,  ECL 

ECL 

PECL 

LVDS/SCI 

Application 

specific 

Target  Cost: 

S/transceiver 

$320  ($  16/Gbps) 

$150 

($  100/Gbps) 

$600  ($66/Gbps) 

Cost  added  to 
next  line 

$/fiber  + 
connectors 

<$200/ 10m 
cable 

<$ 100/Gbps 

The  POINT  passive  polymer  effort  achieved  an  interconnect  density  of  250  lines  per 
board-edge-inch.  Each  line  could  operate  at  up  to  1  gigabits  per  second  (Gbps).  The  total 
(electrical  and  optical)  insertion  loss  achieved  was  -9  dB.  The  waveguide  elements  were 
aligned  passively.  The  fabrication  process  is  MCM  foundry  compatible. 
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3.  The  Computer  Program  BeamPROP 

The  BeamPROP  software  incorporates  advanced  finite-difference  beam  propagation 
techniques  for  simulation,  and  a  modem  graphical  user  interface  for  ease  of  circuit  layout  and 
analysis.  The  benefit  of  good  design  and  modeling  tools  is  well  known  in  the  electronics 
industry,  where  both  device  and  circuit  simulation  programs,  such  as  PISCES  and  SPICE,  have 
been  instrumental  in  advancing  the  availability  and  use  of  integrated  electronic  circuits. 
BeamPROP  brings  this  important  capability  to  the  photonics  area,  and  can  be  an  extremely  useful 
tool  for  research  and  development  groups  in  both  university  and  industrial  environments. 

The  main  program  contains  a  complete  CAD  layout  system  for  design  of  waveguide 
devices  and  circuits,  and  controls  simulation  features  such  as  numerical  parameters,  input  field, 
and  display  and  analysis  options.  The  simulation  program,  which  can  be  executed  within  the 
main  program  or  as  a  standalone,  performs  the  actual  simulation  and  provides  a  graphical  display 
of  the  field  and  other  quantities  of  interest  for  analysis. 

The  CAD  system  was  designed  specifically  for  photonic  circuits,  and  incorporates  a 
natural,  object-oriented  input  model.  Several  unique  features  such  as  user  programmable  circuits 
and  optical  pathways  and  monitors  allow  for  rapid  parameter  variation  and  ease  of  analysis.  The 
CAD  system  allows  a  smooth  transition  to  the  simulation  program,  with  quick  startup  and  easy 
access  to  features.  The  program  provides  intelligent  choices  for  default  numerical  parameters, 
which  are  problem-dependent.  The  simulation  incorporates  integrated,  real-time  analysis  and 
display  of  the  simulation  results,  and  can  be  run  in  the  background.  Both  2D  (one  transverse,  one 
longitudinal)  and  3D  (two  transverse,  one  longitudinal)  simulation  techniques  are  available. 

One  of  the  key  tasks  in  the  POINT  program  was  to  develop  a  software  package  capable  of 
designing  and  analyzing  multimode  waveguide  devices  for  optical  interconnects.  This  outline 
some  of  BeamPROP’s  capabilities  and  applications. 


9 


The  software,  BeamPROP,  was  developed  at  Columbia  University  and  licensed  in  1994 
to  RSoft,  Inc.,  a  commercial  software  firm.  During  the  POINT  program,  Columbia  University 
worked  with  RSoft  to  extend  BeamPROP  so  that  it  would  meet  the  program's  specific  needs  for 
multimode  devices  and  for  interconnects.  The  enhanced  BeamPROP  software,  together  with  the 
multimode  interconnect  device  templates,  forms  the  CAD  tool  developed  under  the  POINT 
program.  The  software  is  capable  of  addressing  such  problems  as: 

Fiber  to  waveguide  coupling 

Source  to  waveguide  and  waveguide  to  detector  coupling 

Waveguide  transitions 

Bend  and  routing  analysis 

Transmission/loss  calculation 

Signal-to-noise,  modal  noise,  and  BER  estimation 

3.1  BeamPROP  Enhancements 

BeamPROP  software  is  a  general  purpose  waveguide  design  package.  At  the  start  of  the 
POINT  program,  the  software  supported  the  following  capabilities: 

CAD  layout,  simulation,  and  design  of  optical  waveguide  devices  and  circuits 
Advanced  parametric  CAD  tool 

Graphical  interface:  input  material  index  parameters,  waveguide  components,  mode 
launch  properties,  etc. 

2D  and  3D  simulation  capability 

Enhanced,  wide-angle  beam  propagation 

Integrated  dynamic  analysis 

Export  of  standard  CAD  files 

PC  or  workstation  based  (Windows  or  UNIX) 


10 


As  the  program  progressed,  Columbia,  working  with  RSoft  made  the  following  additions 
and  improvements  to  the  software  in  response  to  the  POINT  team  requirements: 

Extended  multimode  launch  options 
Enhanced  automatic  variable  scanning  for  random  sampling 
Added  statistical  analysis  tools  for  computing  modal  noise 
Added  bit  error  rate  (BER)  estimation 

Added  a  wide-angle  transparent  boundary  condition  for  multimode  problems 

Added  vector  BPM  capability  for  handling  polarization 

Expanded  its  capability  for  simulating  bends 

Added  new  optimizations,  taking  advantage  of  multimode  structures 

Enhanced  the  wide-angle  BPM  for  multimode  problems 

Improved  the  interface  conditions  improving  the  accuracy  of  high  index  contrast 
Added  new  options  and  utilities  for  modal  decomposition  and  construction  of  input 
fields 

Enhanced  the  batch  capabilities 
Improved  the  file  export  capabilities 

3.2  Theoretical  Background 

The  objective  of  BeamPROP  is  to  provide  a  general  simulation  package  for  computing 
the  propagation  of  light  waves  in  arbitrary  waveguide  geometries.  This  is  a  complex  problem,  in 
general,  and  several  assumptions  are  made  at  the  outset  (many  of  which  are  subsequently 
relaxed).  We  first  limit  our  attention  to  monochromatic  scalar  waves,  which  are  described  by  the 
Helmholtz  equation: 

.,2 

V9  ,  CJL)  9  ,  - 

+  — rn~0  =  0 
c 
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Here,  the  optical  field  (one  component  of  E  or  H)  has  been  written  as  0  (x,y,z)  e~ih  where  0  is 

the  angular  frequency  of  the  light.  The  refractive  index  distribution  in  the  domain  of  interest  is 
assumed  to  be  given  by  n(x,y,z),  and  is  determined  by  the  geometry  of  the  waveguide  circuit. 

The  BeamPROP  software  determines  an  approximate  solution  to  the  above  equation 
using  several  techniques.  The  computational  core  of  the  program  is  based  on  a  finite  difference 
beam  propagation  method  as  described  in  the  paper  by  R.  Scarinozzino  and  R.M.  Osgood,  Jr. 
[75],  and  references  therein.  This  technique  uses  finite  difference  methods  to  solve  the  well- 
known  parabolic  or  paraxial  approximation  of  the  Helmholtz  equation,  which  is  written  as: 

2 ikuz  +  V£u  +  (P  -k2)u=  0 


In  this  equation  we  have  assumed  that  the  wave  is  propagating  primarily  along  z,  and  have  factored  out 
the  rapid  variation  in  the  field  by  writing  <j)(x,y,z)  =  u(x,  y,  z)e'*~ ,  where  k  is  a  constant  representing 
the  characteristic  wavenumber  for  the  propagation.  We  denote  the  wavenumber  in  free  space  by 
k0  =  col c  —  2k / X ,  and  have  introduced  the  notation  k(x,  y,  z)  =  k0n(x,  y,  z)  to  represent  the 

spatially  dependent  wavenumber.  The  above  equation  is  solved  numerically  using  an  implicit 
finite  difference  scheme  as  described  in  reference  1 .  In  addition,  the  program  uses  "transparent" 
boundary  conditions  following  the  work  of  G.R.  Hadley  [76]. 

The  fundamental  physical  limitation  of  the  above  approach  results  from  the  parabolic 
approximation  to  the  Helmholtz  equation,  which  implies  a  paraxiality  condition  on  the  primary 
direction  of  propagation.  These  limitations  can  be  reduced  using  more  accurate  approximations 
to  the  Helniholtz  equation  as  outlined  by  G.R.  Hadley  [77],  BeamPROP  has  the  option  of 
implementing  this  technique,  and  includes  (1,0),  (1,1),  (2,2),  (3,3),  and  (4,4)  Pade7 
approximations. 


BeamPROP  utilizes  vector  beam  propagation  techniques  to  overcome  the  limitations  of  the 
assumption  of  scalar  waves  which  prevents  polarization  effects  from  being  considered.  Those 
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methods  are  based  in  part  on  the  approach  described  by  W.P.  Huang  and  C.L.  Xu  [78, 79]  and 
related  references. 

BeamPROP  uses  a  bidirectional  beam  propagation  (BPM)  algorithim.  The  algorithim 
considers  coupled  forward  and  backward  traveling  waves,  and  can  account  for  reflection 
phenomenon,  including  resonant  effects  as  found  in  grating  structures. 

The  physical  propagation  problem  requires  two  key  pieces  of  information: 

1 .  The  refractive  index  distribution,  n(x,  y,  z)  ■ 

2.  The  input  wave  field,  u(x,  y,  z  =  0) . 

From  these,  the  physics  dictate  the  wave  field  throughout  the  rest  of  the  domain,  u(x,  y,z>  0) . 
The  software  provides  a  way  to  specify  this  information. 

The  solution  algorithm  requires  additional  input  in  the  form  of  numerical  simulation 
parameters  such  as: 

1 .  A  finite  computational  domain  {  x  e  ) } ,  { y  e  ( y^ ,  y^ ) } ,  and 

{Z£  (^min  » ^max  ) ) 

2.  The  transverse  grid  sizes,  Ax  and  Ay. 

3.  The  longitudinal  step  size,  Az. 

The  software  attempts  to  estimate  appropriate  values  for  these  parameters,  but  allows  the  user  to 
override  them.  As  with  any  simulation,  confidence  in  the  accuracy  of  the  numerical  solution 
requires  experimentation  to  determine  the  sensitivity  to  the  numerical  parameters. 
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4. 


BEAMPROP  RESULTS 


Calculated  results  for  passive  waveguides  using  the  CAD  software  BeamPROP 
are  presented  in  this  chapter.  The  work  on  “crosstalk  in  the  optical  backplane”  and 
“simulations  on  multimode  S-bends”  was  conducted  on  the  POINT  program  by  GE  [80], 
The  GE  work  is  based  on  the  use  of  polymers.  One  of  the  polymers  used  on  POINT  is 
Ultem™  which  is  a  trademark  of  the  General  Electric  Company.  Another  polymer  used 
on  the  point  program  was  Polyguide™  which  is  a  trademark  of  Allied  Signal.  Zero-gap, 
single-mode,  directional  couplers  simulations  were  provided  by  the  Air  Force  Research 
Laboratory  [81].  The  single-mode  work  is  based  on  the  use  of  gallium  nitride  (GaN). 

4.1  Calculation  of  crosstalk  in  the  optical  backplane 

During  the  POINT  effort,  GE  performed  a  series  of  simulations  to  determine  an 
upper  bound  on  the  expected  level  of  crosstalk. 

GE  considered  the  basic  problem  of  crosstalk  between  50  pm  wide  waveguides  in 
a  waveguide  array  with  a  100  pm  pitch.  There  are  two  sources  of  crosstalk:  evanescent 
coupling  and  radiation.  In  Figure  2,  the  simulation  results  for  an  array  of  five 
waveguides  using  a  specific  multimode  combination  launched  into  the  center  guide  with 
no  offset  are  shown.  Since  there  is  no  transition  here,  no  radiation  is  expected,  only 
evanescent  coupling.  The  results  indicate  that  the  evanescent  coupling  over  a  1  cm 
length  is  less  than  60  dB,  and  is  therefore  negligible. 

Next,  a  situation  where  the  input  source  is  misaligned  by  an  offset  with  respect  to 
the  center  waveguide  is  investigated.  The  results  for  offsets  of  1,  2,  and  5  pm  (see  Figure 
3)  show  that,  as  the  radiation  induced  by  the  transition  passes  through  the  guide,  the 
crosstalk  level  in  each  guide  reaches  a  maximum  at  a  certain  distance  from  the  input. 
Beyond  that  point,  the  crosstalk  level  decays  slowly.  The  peak  crosstalk  level  ranges 
from  about  -30  dB  for  a  1  pm  offset,  to  about  -20  dB  for  a  5  pm  offset,  giving  one  an 
indication  of  the  worst  case  scenario. 
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Figure  2.  Simulations  for  an  array  of  five  waveguides. 


Figure  3.  Results  for  offsets  of  1 , 2,  and  5  |im. 


With  regard  to  these  results,  there  are  two  important  points  to  note: 


1 .  The  crosstalk  level  is  not  artificial.  Instead,  the  crosstalk  level  is  due  to  light 
trapped  between  the  waveguides,  which  leaks  into  the  adjacent  waveguides. 

2.  These  are  2D  simulations.  One  expects  the  crosstalk  levels  to  be  better  with 
3D  simulations  because  there  is  an  additional  dimension  in  which  to  diffract. 

Simulations  showed  the  beam  diffracts  rapidly  over  the  1  cm  length,  so  we  expect 
that  the  laterally  trapped  light  will  dissipate  considerably  due  to  vertical  diffraction. 
However,  when  the  vertical  field  distribution  was  more  like  the  fundamental  mode, 
diffraction  was  weak  and  would  not  have  been  sufficient  to  dissipate  the  trapped  light. 

4.2  Simulations  On  Multimode  S-Bends 

Beam  propagation  simulations  were  performed  for  2.5°  and  7.5°  S-bends  in 
Polyguide™  to  determine  the  expected  loss  values.  The  simulation  results  for  2.5°are 
shown  in  Figure  4  and  indicate  a  0.85  dB  loss  for  the  full  structure  (two  angular  bends). 
Figure  5  is  a  simulation  that  assumes  7.5°  S-bends  and  indicates  a  2.95  dB  loss  for  the 
full  structure.  For  the  simulations  that  are  represented  by  Figures  4  and  5,  abrupt  S-bends 
(also  called  Z-bends)  were  used.  Note  that  this  is  the  worst  case  scenario. 


Next,  abrupt  S-bends  (also  called  Z-bends)  are  compared  with  that  of  smooth  S- 
bends,  both.  The  1.5°  bend  angle  is  the  largest  angle  that  would  be  needed  in  the  POINT 
receiver  module  where  the  waveguides  would  be  used.  Having  studied  this  S-bend 
theoretically  in  its  abrupt  form,  the  smooth  version  was  modeled.  The  simulation  results 
are  shown  in  Figure  6. 
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Figure  5.  Simulation  results  for  7.5°  S-bends. 


The  smooth  S-bend  simulation  shown  in  Figure  6  reveals  a  bending  loss  of  only 
0.38  dB.  This  confirms  that  we  can  considerably  improve  the  loss  results  by  using  a  non- 
abrupt  design.  Note  that  the  power  loss  in  the  abrupt  structure  occurs  in  two  large  drops 
at  the  bending  junctions,  whereas  the  power  loss  in  the  smooth  structure  occurs  in  two 
small  drops  at  the  first  bend  and  then  at  the  midpoint  (the  point  of  maximum  curvature). 
Thus,  for  a  7.5°  effective  angle  circular  arc,  the  S-bend  angle  is  15°. 


Figure  6.  Simulation  results  for  smooth  7.5°  S-bends. 
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4.3  Simulations  for  single-mode  directional  couplers 

The  POINT  program  utilized  multimode  passive  polymer  waveguides. 
Multimode  waveguides  do  not  have  the  demands  on  lithography  and  processing  that 
single-mode  waveguides.  Passive  polymers  require  individual  lasers  for  every 
waveguide  in  order  to  provide  modulated  information.  Laser  diodes  are  low  impedance 
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current  driven,  not  voltage,  driven  devices.  As  the  modulation  drive  frequency  increases, 
the  demands  on  the  electronic  modulator  become  formidable  because  of  impedance 
matching  requirements  and  required  drive  power.  Optically  active  waveguides  can  act  as 
modulators.  Waveguide  modulators  are  voltage  driven  devices. 

GaN  is  actively  being  researched  as  a  high  temperature  electronic  material  for 
amplifiers.  GaN  is  an  optically  active  material  with  an  E-0  coefficient  for  r^  of  1.9 
pm/V,  for  r33  of  0.5  pm/V,  and  for  rsi  of  3.1  pm/V.  While  it  has  not  been  demonstrated, 
both  electronic  and  optical  devices  may  be  made  on  a  common  substrate  using  the  GaN 
material.  Single-mode  GaN  directional  coupler  switches  were  modeled  using 
BeamPROP  at  the  Air  Force  Research  Laboratory.  For  an  isotropic  medium,  only  the  TE 
mode  is  propagated.  Figure  7  shows  the  power  transfer  between  two  waveguides.  The 
power  can  be  seen  to  transfer  from  the  input  waveguide  to  the  adjacent  waveguide  and 
back  to  the  input  guide  as  a  function  of  distance.  By  controlling  the  total  length  of  the 
waveguide,  one  can  decide  from  which  waveguide  the  optical  energy  will  exit. 
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Figure  7.  Simulation  of  a  single-mode  directional  coupler. 


Pathway, 

Monitor: 

1,  Power’ 

2,  Mode  0 

3,  ModeO 

4,  Neff 

5,  Neff 

6,  Neff 


Active,  single-mode,  waveguide  modulators  are  discussed  in  chapter  5. 


5.0  DESIGN  AND  FABRICATION  OF  NON-LINEAR  OPTICAL  POLYMER 
INTEGRATED  OPTICAL  DEVICES 

The  POINT  program  utilized  multi-mode  passive  waveguides.  The  multimode 
waveguides  could  be  manufactured  with  General  Electric’s  High  Density  Interconnect 
(HDI)  process.  The  HDI  process  does  not  have  the  lithographic  resolution  necessary  to 
realize  effective  single  mode  waveguides.  In  the  chapters  that  follow,  single-mode 
waveguides  will  be  investigated.  Single  mode  waveguides  are  needed  for  active  devices. 

5.1  Electro-Optic  Properties  Of  Nonlinear  Optical  Polymers 

Nonlinear  optical  polymers  in  their  bulk  crystalline  form  typically  behave  like 
biaxial  crystals  [15,  69, 70,  82].  However,  if  one  takes  care  to  ensure  that  the  material  is 
either  amorphous  or  that  the  microcrystalline  dimensions  are  very  small,  undesirable 
propagation  velocity  effects,  which  can  introduce  phase  errors,  will  be  minimized.  The 
NLO  polymer  can  then  be  treated  as  a  uniaxial  crystal  [1 1,  15, 83],  It  has  been  shown  that 
metal  particulates  embedded  in  a  dielectric  material  do  not  show  any  directional  preference 
or  change  in  dielectric  coefficient  in  the  absence  of  an  electric  field,  when  the  particulate 
size  is  much  less  than  the  separation  [83],  Application  of  an  electric  field  in  the  direction  of 
a  chosen  axis  establishes  uniaxial  behavior.  Care  must  be  taken  to  ensure  that  the  material 
does  not  go  through  a  crystalline  material  phase  during  processing.  Uniform  growth  and 
minimized  stress,  as  well  as  maximizing  spontaneous  nucleation  sites,  are  crucial.  Realizing 
a  uniaxial  crystal  structure  using  NLO  polymers  is  material  processing  dependent. 

For  our  analysis,  we  shall  assume  an  NLO  polymer  belonging  to  the  tetragonal 
symmetry  class  4  [15,  82].  Uniaxial  crystals  possess  EO  properties  by  which  a  polarization 
and,  hence,  a  birefringence  can  be  induced  when  an  electric  field  (E)  is  applied  to  the  crystal 
(Pockels  effect)  [26,  29].  The  Pockels  effect  is  a  linear  EO  effect;  the  induced  birefringence 
is  directly  proportional  to  the  first  power  of  the  applied  electric  field.  Thus,  uniaxial  crystals 
exhibit  anisotropic  behavior. 

Figure  8  provides  the  reference  coordinates  for  uniaxial  crystals. 
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Figure  8.  Index  Surfaces  for  (a)  a  Positive  Uniaxial  Crystal  and  (b)  a  Negative  Uniaxial 
Crystal  [29]. 

In  an  anisotropic  crystal,  the  polarization  induced  by  E  is  dependent  on  the  direction 
of  E  [26, 29].  Therefore,  the  speed  at  which  light  travels  through  an  anisotropic  crystal  is 
dependent  upon  the  direction  of  propagation.  A  light  beam  propagating  through  an 
anisotropic  crystal  along  the  optic  axis  (O)  is  defined  as  an  ordinary  ray  [26, 29].  The 
refractive  index  seen  by  this  propagating  beam  is  defined  as  the  ordinary  refractive  index 
(n  y  or  n0).  n  ||  is  independent  of  the  direction  of  the  polarization  [26, 29].  A  light  beam 
propagating  through  an  anisotropic  crystal  at  any  other  direction  to  O  is  defined  as  an 
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extraordinary  ray  [29].  The  refractive  index  seen  by  this  propagating  beam  is  defined  as  the 
extraordinary  refractive  index  (nx  or  ne).  nx  is  dependent  on  the  direction  of  polarization 
[26, 29].  Hence,  nx  is  dependent  on  the  direction  E.  For  our  uniaxial  crystal,  nz  =  n||  and  nx 
=  ny  =  nx.  The  index  surfaces  can,  therefore,  be  represented  by  a  sphere  with  radius  n  \\  an 
ellipsoid  with  radius  n|j  parallel  to  O,  and  radius  nx  normal  to  the  O  [29]  (see  Figure  8).  For 
nx  >  n  || ,  the  uniaxial  crystal  is  known  as  a  positive  type  [29].  For  nx  <  n  || ,  the  uniaxial 
crystal  is  known  as  a  negative  type  [29].  The  sphere  and  ellipse  are  tangent  to  each  other 
only  at  the  points  along  the  optic  axis.  The  amount  of  induced  birefringence  is  defined  as 
the  difference  between  these  two  refractive  indices  and,  as  stated  above,  is  linearly 
proportional  to  E  [26, 29].  Pockels  effect  is  given  mathematically  by  [26, 29] 


n||  -nx  =  mE,  q; 

where  m  is  a  proportionality  constant.  Modulating  the  electric  field  allows  one  to  modulate 
the  input  signal  and  provide  the  basis  for  building  an  electro-optic  switch. 

In  the  next  section  we  will  determine  the  induced  birefringence  and  derive  an 
expression  for  electro-optic  modulation  using  the  EO  properties  of  the  tetragonal  symmetry 
class  4  NLO  polymer. 

5.1.1  Linear  Electro-Optic  Effect 

We  will  begin  our  analysis  with  the  constitutive  equation  [83,  84] 

D  =  £0E  +  P.  (2) 

In  Eq.  (2),  D  is  the  electric  displacement  or  electric  flux  density  vector,  E  is  the  electric  field 
vector,  P  is  the  electric  polarization  vector,  and  £0  is  the  permittivity  of  vacuum.  In  many 
isotropic  materials  the  induced  polarization  is  directly  proportional  to  the  applied  field 
strength,  except  for  the  case  of  very  high  fields.  We  can  write  [82]: 

p  =  £o%E,  (3) 
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where  %  is  the  dielectric  susceptibility  of  the  medium.  Combining  Equations  (2)  and  (3) 
yields  [29,  82] 


D  =  eE, 


(4) 


where  e  =  £o(l  +  %)  is  the  permittivity  of  the  medium. 

The  analysis  is  similar  in  anisotropic  materials.  Equation.  (3)  can  be  expressed  in  its 
tensor  form  as  [15, 82] 


Px  —  £o(XllEx  +  Xl2Ey  XbEz), 

(5) 

Py  =  £o(X2lEx  +  X22Ey  +  fojEJ 

(6) 

Pz  =  £o(X3lEx  +  X32Ey  +  X33Ez), 

(7) 

where  Xij  are  the  coefficients  of  a  3x3  array  called  the  dielectric  susceptibility  tensor  of  the 
medium  for  which  i  =  1, 2  or  3,  and  j  =  1, 2  or  3.  The  suffixes  i  and  j  designate  which  axes 
the  coefficients  correspond  to;  1  represents  x,  2  represents  y  and  3  represents  z.  Choosing  x, 
y  and  z  so  that  the  diagonal  terms  vanish  gives  the  principal  dielectric  axes  of  the  crystal 
[15]: 


Px  —  £o%i  lEx, 


(8) 


Py  —  £<>X22Ey  (9) 

and 

Pz  =  £oX33E2.  (10) 
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Combining  Eq.  (2)  with  Equations  (8)  -  (10)  yields  [15, 29, 82] 


dx  -  £iiEx, 

(11) 

dy  =  822Ey 

and 

(12) 

dz  =  £33^, 

(13) 

where 

£n=£o(l  +%n), 

(14) 

£22  =  £o(l  +X22) 

and 

(15) 

£33  =  £o(l+  X33) 

(16) 

are  the  diagonal  coefficients  of  a  3x3  array  ey, 

of  the  medium. 

which  is  referred  to  as  the  permittivity  tensor 

If  we  introduce  the  phase  velocity  (d  =  (jxe)"1/2),  where  |i  is  the  permeability  of  the 
medium,  it  becomes  evident  that  for  anisotropic  crystals,  u  depends  on  the  direction  of  the 
applied  field.  Thus  we  have  [15,  83,  84] 

'Ox  =  (|J£ii)'1/2, 

(17) 

Dy  =  (jlis22)_1/2 

and 

(18) 

Uz  =  (M£33)'1/2. 

(19) 

For  an  anisotropic  crystal,  at  least  one  of  the  permittivity  coefficients  is  not  equal  to  the 
other  two  (i.e.,  £n  E22  or  833).  Only  the  polarization  parallel  to  the  applied  electric  field 
and,  hence,  the  permittivity  parallel  to  E,  will  contribute  to  the  phase  velocity.  For  example, 
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if  an  input  beam  is  propagating  along  the  y  axis  and  E  is  applied  parallel  to  the  z  axis,  £33 
will  be  the  permittivity,  a  polarization  (Pz)  will  be  induced,  and  the  phase  velocity  will  be  \)z 
(see  Equation  (19)). 

The  index  of  refraction  is  defined  as  the  reciprocal  of  the  relative  phase  velocity  [29, 

83, 84], 

n  =  c/i),  (20) 


where  c  is  the  speed  of  light  in  vacuum. 

Going  to  tensor  form,  we  can  write  the  principal 

refractive  indices  as  [15, 83] 

nx  —  C/\)x  ~  (Pr£rll)  , 

(21) 

ny  =  C/Dy  =  (Pr£r22)1/2 

(22) 

and 

nz  =  c/t)y  =  (Pr£r33)  • 

(23) 

In  Equations  (21)  -  (23),  pr  =  (i/jio  is  the  relative  permeability  of  the  medium,  |io  is  the 
permeability  of  vacuum  and  £nj  =  £i/£0  are  the  coefficients  of  a  3X3  array  called  the  relative 
permittivity  tensor  of  the  medium.  For  nonconductive  media  the  relative  permeability  (pr) 
can  be  assumed  unity  [83],  and  the  principal  refractive  indices  then  become  [15, 83] 


nx  =  (£ii/£o)1/2, 

(24) 

ny  =  (£22/£o)1/2 

(25) 

nz  =  (£33/e0)1/2- 

(26) 
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Using  the  definition  for  the  refractive  index  along  with  the  analysis  for  anisotropic  crystals 
given  above,  it  becomes  evident  that  one  can  induce  a  change  in  the  refractive  index  with  an 
applied  electric  field.  This  change  will,  in  turn,  induce  birefringence. 

5.1.2  Index  Ellipsoid 

We  will  start  with  the  general  index  ellipsoid  for  a  crystal  in  its  simplest  form.  This 
will  allow  us  to  determine  the  direction  of  polarization  as  well  as  the  corresponding 
refractive  indices  of  the  crystal.  The  index  ellipsoid  is  given  in  the  principal  coordinate 
system  as  [27, 29, 85] 

x2/nx2  +  y2/ny2  +  z2/nz2  =  1 .  (27) 

In  Equation  (27),  x,  y  and  z  are  the  principal  axes  (the  direction  in  the  crystal  along  which 
the  electric  field  vector  (E)  and  the  electric  displacement  vector  (d)  are  parallel).  l/nx2,  l/ny2 
and  l/nz"  are  the  principal  refractive  indices  along  their  respective  axes.  The  lengths  of  the 
major  axes  of  the  index  ellipsoid,  represented  by  Equation  (27),  are  2nx,  2ny  and  2nz,  parallel 
to  directions  x,  y  and  z,  respectively  [15,  85].  The  index  ellipsoid  is  used  to  find  the  two 
refractive  indices  and  two  corresponding  directions  of  d  associated  with  the  two 
independent  plane  waves  propagating  along  an  arbitrary  direction  s  in  the  crystal.  These 
indices  are  determined  using  Equation  (27).  Visually  depicted  in  Figure  9,  the  major  and 
minor  axes  of  the  ellipsoid  are  in  the  directions  dj  and  d2,  respectively,  with  the  length  of  di 
equal  to  2m,  and 


Figure  9.  Inner  Ellipse  is  the  Intersection  of  Index  Ellipsoid  with  the  Plane  Perpendicular  to 
s  [15]. 
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Next  we  will  introduce  the  optical  dielectric  impermeability  tensor  of  a  crystal.  The 
coefficients  (riy)  of  this  tensor  depend  on  the  distribution  of  bond  charges  in  the  material 
[15, 82],  riy  are  found  by  taking  the  reciprocal  of  the  relative  permittivity  or  dielectric 
constant  [82].  riy  has  been  defined  in  terms  of  the  refractive  index  of  the  crystal  as  [82] 

tlij =  £o/£ij =  1/ny  (28) 

(see  Equations  (24)  -  (26)).  Application  of  an  electric  field  E  to  the  crystal  results  in  a 
redistribution  of  the  bond  charges  of  the  crystal,  and  a  possible  deformation  of  the  ion 
lattice.  This  causes  a  change  in  the  optical  impermeability  tensor  and,  hence,  a  change  in  the 
refractive  index  of  the  crystal.  This  allows  us  to  write  the  index  ellipsoid  in  the  presence  of 
an  electric  field  as  [15]: 

(l/nx2)x2  +  (l/ny2)y2  +  (l/nz2)z2  +  Ajx2  +A2y2  +  A3Z2  +  A4yz  +  A5ZX 
+  A6xy  +  Kjzy  +  Asxz  +  Agyx  =  1 .  (29) 

For  E  =  0,  A]  through  A9  =  0,  and  Equation  (29)  reduces  to  Equation  (27).  Introducing  the 
notations  for  subscripts: 

x2  =  xx  =  (ll)->  1 
y2  =  yy  =  (22)^2 
z2  =  zz  =  (33)  — »  3 
yz  =  (23)  ->  4 

zx  =  (31)  — »  5  (30) 

xy  =  (12)  — >  6 
zy  =  (32)  -»  7 
xz  =  (13)  -»  8 
yx  =  (21)  — >  9 
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And, 


(l/n2)i  =  l/nx2  =  Aj 
(l/n2)2=l/ny2  =  A2 
(l/n2)3=l/nz2  =  A3 
(l/n2)4  =  l/nyn2  =  A4 

(l/n2)5=  l/nznx  =  A5  (31) 

(1/n  )6 =  l/nxny  —  Ag 

(1/n  )7  =  l/nzny  =  A7 

(l/n2)8=  l/nxnz  =  A8 

(l/n2)9  =  l/nynx  =  A9, 


Equation  (29)  becomes 

[l/nx2  +  (l/n2)i]x2  +  [l/ny2  +  (l/n^Jy2  +  [l/nz2  +  (l/n2)3]z2 
+  [(l/n2)4]yz  +  [(l/n2)5]zx  +  [(l/n2)6]xy  +  [(l/n2)7]zy 

+  [(l/n2)8]xz  +  [(l/n2)9]yx  =  1.  (32) 

It  has  been  shown  that  riy  is  symmetric,  provided  that  the  medium  is  lossless  and 
optically  active  [15].  This  allows  the  permutation  of  i  and  j,  and  the  notation  presented  in 
Equations  (31)  and  (32)  are  reduced  to: 

x2  =  xx  =  (ll)-»  1 

y2  =  yy  =  (22)  2 

z2  =  zz  =  (33)  -4  3  (33) 

yz  =  zy  =  (23)  =  (32)  — »  4 
zx  =  xz  =  (31)  =  (13)  — >  5 
xy  =  yx  =  (12)  =  (21)->6. 
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with: 


(1/n2),  =  l/nx2 
(l/n2)2  =  l/ny2 

(l/n2)3  =  l/nz2  (34) 

(l/n2)4  =  l/nynz  =  l/nzny 
(l/n2)5  =  l/nznx  =  l/nxnz 
(l/n2)6  =  l/nxny  =  l/nynx. 

Substituting  Equations  (33)  and  (34)  into  Equation  (32)  yields  [27] 

[l/nx2  +  (l/n2)i]x2  +  [l/ny2  +  ( l/n2)2]y2  +  [l/nz2  +  (l/n2)3]z2 

+  2[(l/n2)4]yz  +  2[(l/n2)5]xz  +  2[(l/n2)6]xy  =  1 .  (35) 

5.1.3  Electro-Optic  Coefficients 

The  electro-optic  (EO)  coefficients  are  traditionally  defined  as  [15,  29] 

fiij(E)  ■  Pii(O)  AT|jj  =  tykEfc  +  SykjE^E].  (36) 

In  Equation  (36),  E  is  the  applied  electric  field,  are  the  coefficients  of  a  3x9  array  called 
the  linear  electro-optic  (Pockels)  tensor  of  the  medium.  Syki  are  the  coefficients  of  a  6x6 
array  called  the  quadratic  electro-optic  (Kerr)  tensor  of  the  medium,  for  which  k  =  1, 2,  or  3 
and  1  =  1, 2,  or  3.  In  the  equation  above,  we  have  omitted  terms  higher  than  the  quadratic. 
These  higher  order  effects  are  so  weak  compared  with  the  linear  and  quadratic  effects,  that 
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they  can  be  neglected  for  most  applications  [15,  29].  For  centrosymmetric  crystals  (crystals 
possessing  inversion  symmetry  for  all  physical  properties),  the  linear  EO  effect  vanishes 
[71]  (i.e.,  %  =  0).  This  is  caused  by  the  spatial  inversion  symmetry  of  such  crystals. 
Therefore,  the  quadratic  EO  effect  will  dominate  [15].  For  noncentrosymmetric  crystals 
(crystals  that  do  not  possess  inversion  symmetry  for  all  physical  properties),  both  the  linear 
and  quadratic  effects  will  be  present  [15].  From  the  section  above,  the  EO  effects  will 
depend  on  the  ratio  of  the  applied  electric  field  (E)  to  the  intra-atomic  electric  field  (Eu) 
binding  the  charged  particles.  This  ratio  can  be  estimated  by  finding  the  magnitude  of  Eu  at 
a  position  (T)  arising  from  a  point  charge  (q)  located  at  a  distance  (u)  from  T.  Eu  at  T  is 
given  by  [83,  84] 

Eu  =  (l/47i£o)q/u2.  (37) 

The  minimum  field  inside  the  atom  Eu  can  be  determined  by  substituting  the  lattice  constant 
ao  of  the  crystal  for  u,  and  the  elementary  electron  charge  for  q.  If  the  minimum  electric 
field  inside  the  atom  is  much  larger  than  the  applied  field,  then  the  quadratic  EO  effect 
should  be  small  enough  when  compared  to  the  linear  EO  effect.  This  will  allow  us  to 
neglect  the  Syu  terms.  Equation  (36)  becomes 


fiij(E)  -  T)ij(O)  AT|ij  =  rUkEk. 


(38) 


The  permutation  property,  along  with  Equations  (28),  (33),  (34)  and  (38),  yields  [26, 29] 

3 

A(l/n2)i=^  rikEk,  (39) 

i=l 

where,  for  example,  A(l/n2)5  =  r51Ei  +  r52E2  +  r53E3.  Substitution  of  Eq.  (39)  into  Eq.  (35) 
gives  the  general  index  ellipsoid  in  the  presence  of  an  electric  field  in  terms  of  nk  [15]: 

[l/nx2  +  r[kEk]x2  +  [l/ny2  +  r2kEk]y2  +  [j  /n22  +  r3kEk]z2 

+  2[r4kEk]yz  +  2[r5kEk]xz  +  2[r6kEk]xy  =  1  (40) 
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In  order  to  determine  the  EO  tensor  for  crystals  belonging  to  the  tetragonal  4  class, 
we  will  begin  with  the  EO  tensor  for  the  triclinic  1  crystal  group  in  matrix  form  [27]: 


rn 

ri2 

rn 

r2i 

T22 

r23 

T31 

T32 

T33 

V4\ 

r42 

T43 

rsi 

T52 

T53 

T61 

V62 

T63_ 

Neglecting  the  permutation  property  for  the  moment,  Equation  (41)  expands  to 


(41) 


rill 

rm 

rn3 

r22i 

r222 

r223 

1*331 

T332 

T333 

r23i 

r232 

r233 

T321 

1*322 

T323 

ri31 

ri32 

ri33 

T311 

r3i2 

T313 

rm 

n22 

1*123 

r2n 

r2i2 

r213_ 

(42) 


Since  tetragonal  crystals  have  axes  of  fourfold  symmetry,  we  can  let  z  be  the  axis  of  two 
fold  rotation,  which  can  be  represented  by  the  matrix  [15] 


an  = 


-1 

0 

0 


0 

-1 

0 


0 

0 

1 


(43) 


Equation  (43)  dictates  that  the  suffixes  of  the  EO  coefficients  change  as  follows:  1  — >  -1 , 
2  — >  -2  and  3  — »  3  [82].  Substituting  this  property  into  Equation  (42)  yields 
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where 


Tin 

1*113 

"  V222 

T223 

“  T332 

1*333 

T232 

“T233 

T322 

~  T323 

ri32 

”  ri33 

T312 

T313 

“  ri22 

n?3 

“  T212 

T213_ 

(44) 


r'jk  rni  r(-l)(-l)(-l)  rill,  rijk  —  r123  r(.i)(.2)(3)  —  r123;  rijk  —  I223  — >  r(.2)(-2)(3)  =  1*223 j 

etc. 

Setting  rijk  =  rijk  leads  to  the  following  condition:  rp  =  -r^k,  only  if  rp  =  0  [15].  Equating 
Equations  (42)  and  (44)  with  this  condition  yields 

0  0  no 

0  0  f223 

0  0  1*333 

T231  f232  0 

rijk  =  T321  T322  0  (45) 

ri3i  ri32  0 

T311  T312  0 

0  0  n  23 

.  0  0  T213. 

Since  tetragonal  4  class  crystals  also  includes  a  diad  axis  parallel  to  z,  the  axes  transform  as 
follows  [71]:  1  — »  -2, 2  1  and  3  — >  3.  Hence, 
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where 


0  rn3 

0  n  13 

0  r333 

fl3I  0 

f3Il  ® 

T231  0 

“  T321  0 

0  0 

0  0. 


(46) 


%  =  rj  13  — »  r(.2X-2X3)  —  1*223;  %  =  1*231  r(l)(3X-2)  =  -1*132;  %  =  1*321  — : >  l*(3Xl)(-2)  =  _r312; 

I*ijk  =  1*131  — >  l*(-2X3X-2)  =  1*232;  %  =  1*311  I*(3X-2)-21)  =  1*322;  %  =  1*231  ~ »  r(3Xl)(2)  =  r312; 

%  =  ri23  — » i*(-2xiK3) =  -1*213  and  r,jk  =  r2i3  — >  i*(ix-2)(3>  =  -ri23-  For  ri23  =  -1*213  and  1*213  =  -1*123; 
which  proves  that 

1*123  =  1*213  =  0.  (47) 

Going  back  to  contracted  notation,  the  electro-optic  tensor  (rik)  for  crystals  belonging  to  the 
tetragonal  4  group  can  be  represented  by  the  well  known  matrix  [15] 

0  0  ri3 

0  0  ri3 

0  0  r33 

T41  T51  0 

rsi  -  r  4 1  0 

0  0  0 
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Substituting  the  EO  tensor  for  tetragonal  crystals  back  into  the  index  ellipsoid  (Equation 
(40))  gives  [15] 

[l/nx2  +  ri3Ez]x2  +  [l/ny2  +  r^EJy2  +  [l/n22  +  r33Ez]z2 

+  2yzr41Ex  +  2yzr5iEy  +  2xzr5iEx  -  2xzr4iEy  =  1 .  (49) 

Taking  into  account  that  tetragonal  4  crystals  are  uniaxial  allows  us  to  write  [85] 

[l/n02  +  ri3Ez]x2  +  [l/n02  +  r^EJy2  +  [l/ne2  +  r33Ez]z2 

+  2yzr4]Ex  +  2yzr51Ey  +  2xzr5]Ex  -  2xzr4iEy  =  1,  (50) 

where  nD  and  ne  are  the  ordinary  and  exrtraordinary  refractive  indices,  respectively. 

5.1.4  Induced  Refractive  Indices  and  Birefringence 

We  need  to  determine  the  direction  and  magnitude  of  the  principal  axes  of  the  new 
ellipsoid  represented  by  Eq.  (50),  since  they  may  not  coincide  with  x,  y,  z.  For  this  analysis, 
we  are  fabricating  a  transverse  modulated  EO  switch  [15].  This  dictates  that  the  electrodes 
be  placed  parallel  to  the  waveguides  along  the  optic  axis.  Propagation  of  the  input  beam 
will  be  along  the  y  axis  through  the  electro-optic  device  of  length  L.  The  electrodes  will  be 
placed  on  the  top  and  bottom  of  the  coupler  running  along  y  and  separated  by  a  distance  d. 
The  electric  field  E  will  be  applied  along  the  z  axis,  normal  to  y.  Figure  9  illustrates  a 


E 


Figure  10.  Transverse  Electro-Optic  Modulator  using  a  Tetragonal  4  Group  Crystal. 
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transverse  modulator.  Since  the  electric  field  is  in  the  z  direction,  we  have 


Ex  =  Ey  =  0 
and 

EZ  =  E.  C 

Substituting  Equation  (51)  into  Equation  (50)  yields  the  index  ellipsoid 

[l/n02  +  r13E]x2  +  [l/n02  +  r13E]y2  +  [l/ne2  +  r33E]z2.  C 

Since  no  mixed  terms  appear  in  Equation  (52),  the  principal  axes  of  the  new  index 
ellipsoid  remain  unchanged.  Using  Equation  (39),  we  have 


[(l/n02)  +  A(l/n02)]x2  +  [(l/n02)  +  A(l/n02)]y2  +  [(l/ne2)  +  A(l/ne2)]z2  =  1.  (! 

From  Equation  (53)  we  can  write 

l/(nx)2  =  l/n02  +  A(l/n02),  0 

l/(ny)2  =  l/n02  +  A(l/n02)  C 

and 

l/(nz)2  =  l/ne2  +  A(l/ne2).  C 


Using  Equation  (1)  for  Pockels  effect  along  with  Equations  (52),  (38)  and  (39)  yields  the 
induced  change  in  the  refractive  index 

An  =  m'r^E  =  m'A(l/n2),  C 


where  m"  is  a  proportionality  factor  we  will  use  to  equate  An  and  A(l/n  ).  If  we  assume 
rijkE  «  1/n2,  (i.e.,  A(l/n2)  sufficiently  small),  it  allows  us  to  use  the  relation  [15] 


(58) 


A(l/n2)  =  3  (1/n2)  =  -(2/n3)  9  n. 


Multiplying  both  sides  by  -(n3/2)  yields 


3n  =  -(n3/2)9  (1/n2). 


(59) 


Equating  Equations  (57)  and  (59)  leads  to  [15] 


An  =  m’  ~  -(n3/2)A(l/n2) 


(60) 


Using  Equation  (39)  once  again,  we  can  write 


An  =  -(n3/2)rikE. 


(61) 


This  leads  to  the  indices  [15] 

nx  =  n0  +  An0  =  n0  -  (n03/2)r13E, 

(62) 

ny  =  n0  +  An0  =  n0  -  (n03/2)r13E 

(63) 

and 

nz  =  ne  +  Ane  =  ne  -  (ne3/2)r33E 

(64) 

The  induced  birefringence  for  the  new  axes  with  the  input  beam  propagating  perpendicular 
to  the  applied  electric  field  (i.e.,  for  transverse  modulation)  can  now  be  given  by  [15,  26] 

nz  -  nx  =  (ne  -  n0)  -  [(ne3/2)r33  -  (n03/2)r13]E.  (65) 

Equation  (65)  is  the  definition  for  the  induced  birefringence  that  we  will  use  to  describe  the 
phase  retardation  experienced  by  the  input  beam  as  it  propagates  through  the  crystal. 
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5.2  ELECTRO-OPTIC  MODULATORS 


For  our  analysis  we  will  consider  an  optical  field  that  is  incident  and  normal  to  the 
xz  plane  propagating  along  the  y  axis  with  E  parallel  to  the  z  axis.  At  the  input  plane  y  =  0, 
the  optical  field  can  be  resolved  into  two  mutually  orthogonal  components  polarized  along  x 
and  z.  The  x  component  propagates  as  [27, 29] 

Ex  =  Aexp[i(cot  -  (oVc)nxz)],  (66) 

where  co  is  the  angular  frequency,  t  is  time  and  A  is  a  constant.  The  z  component  propagates 
as  [27,  29] 

Ez  =  Aexp[i(cot  -  (co/c)nzy)].  (67) 

The  difference  in  phase  between  these  two  components  at  the  output  plane  y  =  L  is  called 
the  phase  retardation  [27, 29].  It  is  given  by  the  difference  of  the  exponential  terms  in 
Equations  (66)  and  (67).  Performing  the  subtraction  yields 


Txz  =  (oVc)(nz  -  nx)L,  (68) 

where  [nz  -  nx]  is  the  induced  birefringence  found  in  the  previous  section.  Substituting  for 
the  induced  birefringence  using  Equation  (65)  gives  [27, 29] 


rxz  =  (fflfcXGle  -  no)  -  ((ne3/2)r33  -  (n03/2)r13)E)L. 


(69) 


Using  the  relation  co/c  =  2ro9i  with  A  defined  as  the  wavelength,  the  phase  retardation 
becomes  [29, 74] 


rxz  =  (27t/A)(ne  -  n0)L  -  (7t/A)(ne3r33  -  n03r13)VL/d, 


(70) 
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where  V  is  the  voltage  applied  between  electrodes  separated  by  a  distance  d.  Due  to  the 
natural  birefringence  terrn  in  Equation  (70),  a  Babinet  or  Soleil  phase  compensator  will  be 
required.  This  will  allow  us  to  adjust  the  phase  retardation,  that  is  present  in  the  absence  of 
the  applied  voltage,  to  equal  an  odd  integral  multiple  of  tc/2.  We  then  want  a  voltage 
dependant  phase  retardation  of  180°,  or  Txz  =  7i.  The  voltage  necessary  to  realize  a  it  phase 
change  is  defined  as  half  wave  voltage  [27, 29].  Equating  the  voltage  dependant  phase  term 
of  Equation  (70)  to  it  gives  us  [15] 


V*  =  dA/L(ne3r33  -  n03ri3). 


(71) 


Vjc  is  the  voltage  necessary  to  realize  half  wave  polarization  in  tetragonal  crystals. 
Figure  1 1  illustrates  a  Mach  Zehnder  (MZ)  interferometer. 


Figure  1 1 .  Operation  of  a  Mach  Zehnder  Interferometer. 
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interferometer  is  a  form  of  transverse  electro-optic  modulator.  The  input  signal  is  first  split 
equally  into  two  separate  channels.  These  two  channels  are  separated  by  enough  distance  to 
prevent  any  crosstalk  or  coupling.  The  two  channels  are  later  recombined  back  into  one 
channel.  Voltage  is  applied  to  only  one  of  the  separated  channels.  When  voltage  is  absent, 
half  of  the  signal  propagating  in  channel  one  recombines  constructively  with  half  of  the 
signal  propagating  in  channel  two,  and  the  output  equals  the  that  present  at  the  input. 
Voltage  applied  to  one  of  the  channels  results  in  a  phase  shift  of  the  propagating  beam  only 
in  that  channel.  When  enough  voltage  is  applied  to  shift  the  phase  of  the  beam  in  that 
channel  by  tc,  then  half  of  the  signal  propagating  in  channel  one  recombines  destructively 
with  half  of  the  signal  propagating  in  channel  two,  and  the  output  becomes  zero.  By 
periodically  switching  the  voltage  on  and  off,  one  can  realize  sinusoidal  modulation  of  the 
optical  signal. 
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53  DIRECTIONAL  COUPLERS 


The  directional  coupler  switch  is  a  device  consisting  of  parallel  channel  waveguides 
separated  by  a  finite  distance.  The  coupling  between  the  modes  of  the  parallel  waveguides 
results  in  an  exchange  of  power  between  guided  modes  of  adjacent  waveguides  [15].  This  is 
referred  to  as  directional  coupling.  The  eigenmodes  of  the  coupling  region  consist  of  one 
symmetric  and  one  antisymmetric  mode  [15, 24, 25, 27, 28, 86-89].  Treatment  of 
waveguide  coupling  can  be  performed  by  coupled-mode  theory  [15, 27, 28, 86], 

A  of  directional  coupler  is  illustrated  in  Figure  12.  In  the  figure,  s  is  the  separation 
between  the  waveguides. 


) 

Optical 

Input 


)  j 

Symmetric  Antisymmetric 

Mode  Mode 


Figure  12.  Operation  of  a  Directional  Coupler 
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Consider  Pij  as  the  input  power  to  waveguide  1 ,  Pi2  as  the  input  power  to  waveguide 
2,  Po,  as  the  output  power  from  waveguide  1  and  P02  as  the  output  power  from  waveguide 
2.  Let  waveguide  1  propagate  through  the  symmetric  mode  and  waveguide  2  propagate 
through  the  antisymmetric  mode.  Assume  that  the  waveguides  are  not  too  close,  so  that  the 
overlap  integral  of  the  mode  functions  will  be  small  [15, 27,  28,  86].  With  a  single  input  to 
waveguide  1 ,  Pii  =  1  and  Pi2  =  0,  the  output  power  from  waveguides  1  and  2  are  given  as 
[15,  27, 28,  86,  87] 

Poj(0  =  Pii  -  Po 2(/)  =  Pi,(l  -  (i^/Ck2  +  82))(1  -  cos2[(k2  +  82)1/2(/)]))  (72) 

and 

Po2(0  =  PiitOcV  +  82))(sin2[(K2  +  82)1/2(/)])),  (73) 

where  28  =  pi  -  P2  is  the  difference  in  the  propagation  constants  between  two  adjacent 
uncoupled  waveguides  1  and  2,  k  is  the  coupling  constant  between  the  two  adjacent 
waveguides.  /  is  the  interaction  length  over  which  the  two  waveguides  are  coupled.  By 
definition,  k  =  7t/2L,  where  /  =  L  is  the  length  required  for  complete  coupling  of  the  power 
from  one  waveguide  to  the  other  [15-19,  24, 25, 27,  28,  86-90]. 

For  symmetric  waveguides,  the  phase  velocities  in  the  two  modes  are  equal  (i.e.,  pi 
=  P2  or  8  =  0).  For  complete  coupling,  Eqs.  1 .4. 1  and  1 .4.2  become  [  1 5, 27, 28,  86,  87] 


Po,(L)  =  Pi,(l  -  (1  -  cos2[tc/2]))  =  0 

(74) 

Po2(L)  =  Pii(sinV2])  =  Pii. 

(75) 

From  Eqs.  (72)  and  (73),  one  can  see  that  the  input  power  Pii  into  waveguide  1  will  exit 
waveguide  2  as  P02  at  a  distance  L.  Using  Eqs.  (74)  and  (75),  complete  coupling  occurs 
when  L  =  7t/2K. 

Applying  an  electric  field  to  either  waveguide  1  or  waveguide  2  over  the  distance  L„ 
where  the  waveguides  are  coupled,  will  change  that  particular  waveguide's  propagation 
constant  p.  Total  switching  occurs  when  8  =  31/2tc/2L  [15,  86,  87].  This  can  be  verified  by 
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the  substitution  of  5  =  31/27t/2L  into  Eqs.  (72)  and  (73).  By  doing  this,  we  get  [15, 27, 28, 
86, 87] 


Poi(L)  =  Pi,(l  -  (1/4)(1  -  cos2[7t]))  =  Pij 
and 

Po2(L)  =  Pi1((l/4)(sin2[7i]))  =  0. 

Using  the  definitions  for  8  and  k  given  above,  one  gets  total  switching  to  occur  when  [15, 
27,28, 86] 

(Pi  -  P2)  =  31/27t/L.  (' 

If  we  apply  an  electric  field  (E)  over  waveguide  2,  for  instance,  this  will  induce  a 
change  in  the  refractive  index  of  waveguide  2  by  an  amount  proportional  to  E. 
Mathematically,  this  is  given  by  [15,  27-29, 86] 


(76) 

(77) 


An2  =  (n23/2)rE,  (79) 

where  n2  is  the  refractive  index  of  waveguide  2  and  r  is  the  electro-optic  coefficient  of  the 
waveguiding  material.  To  apply  an  electric  field  to  one  of  the  waveguides,  we  place 
electrodes  on  the  top  and  bottom  of  the  waveguide  over  the  coupling  length  (L)  and  apply  a 
voltage  (V)  to  the  top  electrode.  The  electric  field  is  then  related  to  the  applied  voltage  by  V 
=  Ed,  where  d  is  the  distance  between  the  top  and  bottom  electrode.  The  total  index  (n2')  of 
the  waveguide  2  region  between  the  electrodes  can  be  defined  by 

n2'  =  n2  +  An2  =  n2  +  (n23/2)rV/d.  (80) 

By  definition,  P  =  27mA,,  where  X  is  the  wavelength  of  the  source  [15, 27,  28,  86].  If  we  let 
waveguide  1  and  waveguide  2  be  identical  (n!  =  n2  =  n),  and  substitute  Eq.  (80)  into  (78), 
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along  with  the  definition  for  (3,  this  yields  [15, 27, 28,  86] 

(27t/A.)(n'  -  n)(n3/2)(rV/d)  =  (27t/A)(n3/2)(rV/d)  =  3I/27t/L.  (81) 

We  can  now  find  the  voltage  required  to  switch  the  coupling  completely.  By  rewriting  Eq. 
(81),  we  get 

V  =  (31/2Xd)/(n23rL).  (82) 

The  EO  coefficient  r  =  r33  can  be  substituted  for  the  assumed  crystal  orientation  presented  in 
the  previous  section. 

5.3.1  Zero-Gap  Directional  Couplers 

As  with  directional  couplers,  the  interaction  region  of  a  zero-gap  directional  coupler 
will  support  two  guided  symmetric  modes,  and  two  guided  antisymmetric  modes  [14, 32, 
91-96].  The  difference  between  them,  however,  results  from  the  fact  that  we  have  zero 
spacing  between  the  two  symmetric  waveguides  in  a  zero-gap  directional  coupler.  These 
modes  will  now  be  propagating  in  a  single  waveguide.  We  will  again  have  a  voltage 
independent  phase  difference  that  accumulates  over  the  interaction  length  (the  two  mode 
section),  due  to  the  difference  in  the  phase  velocities  of  the  two  orthogonal  modes,  and  a 
phase  difference  over  the  two  mode  section,  brought  about  by  the  voltage  induced  change  in 
refractive  index  (EO  effect)  [15, 26,  82-84]. 

The  electro-optic  dependent  phase  difference  will  be  denoted  by  (ps  -  Pas)/2  =  Ap/2, 
where  ps  and  P^  are  the  propagation  constants  of  the  symmetric  and  antisymmetric  modes, 
respectively  [92].  The  value  of  AP  required  for  complete  coupling  of  energy  from  the  top 
half  of  the  two  mode  waveguide  to  the  bottom  is  given  by  [92] 

Ap  =  7t/L.  (83) 
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A|3  is  varied  by  changing  the  index  of  refraction  of  the  waveguiding  material  in  the  two¬ 
mode  section.  We  will  denote  the  refractive  index  change,  as  we  did  in  the  last  section,  as 
An.  The  change  in  the  difference  between  the  propagation  constants  of  the  two  modes  can 
then  be  written  as  [92] 

A(Ap)  =  And  (A(3)/9  n.  (84) 

Equating  Eqs.  (83)  and  (84)  yields  [92] 

An  =  (7t/L)/(  3  (Ap)/  3  n).  (85) 

In  the  last  section  it  was  shown  that  for  linear  EO  materials,  the  refractive  index  changes 
proportionally  to  the  applied  voltage.  With  this  relation,  it  can  be  seen  from  Eq.  (85)  that 
a  (AP)/ 3  n  is  inversely  proportional  to  the  applied  voltage,  d  (Ap)/a  n  is  also  strongly 
dependent  on  the  waveguide  separation  s  for  the  directional  coupler  [92].  Papuchon  and 
Roy  [92]  showed  that  as  s  decreases,  d  (Ap)/a  n  increases,  and  as  the  value  of  d  (Ap)/a  n 
increases,  the  voltage  required  for  switching  decreases.  Therefore,  the  maximum  value  of 
a  (Ap)/a  n  and,  thus,  the  minimum  voltage  required  to  switch  the  coupling  guides,  occurs 
when  s  =  0.  When  s  =  0,  we  have  a  zero-gap  directional  coupler  [92],  Experiments  have 
shown  that  the  value  d  (Afi)/d  n  varies  by  an  order  of  magnitude  between  the  conventional 
directional  coupler  and  the  zero-gap  directional  coupler  [92].  This  verifies  that  zero-gap 
directional  couplers  require  less  voltage  to  switch  the  coupling  than  do  conventional 
directional  couplers. 

Expressions  for  the  power  coupling  in  zero-gap  directional  couplers  are  similar  to 
those  used  for  conventional  directional  couplers.  They  are  given  by  [97-99] 

Poi(/)  =  PiiCos2[(Ap/2  +  0>o)/]  (86) 

and 

Po2(/)  =  Piisin2[(Ap/2  +  Oo)/],  (87) 
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where  <E>o  is  the  accumulated,  voltage  independent  phase  difference  in  the  phase  velocities 
of  the  symmetric  and  antisymmetric  modes.  With  no  voltage  applied  (i.e.,  A0  =  0),  coupling 
occurs  when  d>o  =  n/2L.  When  enough  voltage  is  applied  so  that  A(3  =  rc/L,  then  complete 
switching  occurs.  Equating  Ap  with  Eq.  (80)  yields 

(27t/X)(n3/2)(rV/d)  =  nl L.  (88) 

The  voltage  required  for  complete  switching  is  then  given  by 

V  =  (A.d)/(n3rL),  (89) 

where  n  is  the  refractive  index  of  the  interaction  region.  Like  the  directional  coupler,  the  EO 
coefficient  r  =  ^3.  By  comparing  Eq.  (89)  with  Eq.  (82),  one  can  see  that  the  voltage 
required  to  completely  switch  the  optical  input  from  one  channel  to  the  other  is  less  for  the 
zero-gap  directional  coupler.  The  zero-gap  allows  us  to  treat  the  center  interaction  region  as 
one  waveguide.  Operation  of  the  zero-gap  directional  coupler  switch  is  illustrated  in  Fig. 

13.  Another  advantage  of  zero-gap  directional  couplers,  in  addition  to  operating  with  less 
voltage,  is  that  the  interaction  length  is  much  shorter  than  the  interaction  length  of 
conventional  directional  couplers  [24, 32-35,  37,  89,  90].  This  results  in  the  reduction  of  real 
estate  required. 
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Figure  13.  Operation  of  a  Zero-Gap  Directional  Coupler. 


5.3.2  Wavefront  Tilt 


An  additional  parameter  which  needs  to  be  taken  into  account  for  directional 
couplers  and  zero-gap  directional  couplers  is  wavefront  tilt,  and  its  affect  on  waveguide 
coupling  in  the  branching  section.  The  larger  the  tilt  angle  a  of  the  branching  section,  the 
less  the  waveguide  coupling  [91].  Tilt  introduces  phase  differences  that  decrease  waveguide 
interaction  in  the  branching  section  of  waveguide  coupling.  The  maximum  allowable  tilt 
angle  can  be  obtained  by  forcing  the  branching  region  of  waveguide  coupling  to  be  within  a 
single  lobe  of  the  standing  wave  generated  by  the  two  waves  propagating  along  the  tilted 
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waveguide  portions.  This  gives  the  limit  [91] 


sin(a/2)  <  A74nw,  (90) 

where  w  is  the  width  of  the  waveguides. 


5.3.3  Interaction  Length  Of  The  Directional  Coupler 

In  this  section,  the  interaction  length  of  a  directional  coupler  is  determined  by 
substituting  the  refractive  index  of  the  NLO  polymer  into  several  dielectric  strip  waveguide 
directional  coupler  models.  Driving  the  separation  between  the  waveguides  to  zero  allows 
us  to  determine  the  interaction  length  of  the  zero-gap  directional  coupler. 

Directional  coupling  is  the  exchange  of  power  between  the  guided  modes  in  adjacent 
parallel  waveguides  separated  by  a  finite  distance.  Fig.  14  is  a  schematic  of  a  conventional 
directional  coupler  device.  The  index  of  refraction  of  the  core  region  is  n  while  that  of  the 
clad  regions  is  n(  1  -  A).  This  is  less  than  the  index  of  the  core  region  by  an  amount  A.  The 
dimensions  a,  b  and  c  are,  respectively,  the  width,  thickness  and  separation  of  the 
waveguides.  The  interaction  length,  the  distance  required  for  total  exchange  of  power  from 
one  guide  to  the  other,  is  L  [77].  The  wave  is  assumed  to  propagate  along  the  y  axis. 


Figure  14.  Rectangular  Waveguide  Directional  Coupler. 
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Extensive  treatment  of  optical  waveguide  and  directional  coupler  theory  has  been 
addressed  by  many  researchers  for  both  continuous  two-dimensional  film  waveguides  and 
dielectric  strip  waveguides.  Various  methods,  such  as  dielectric  waveguide,  coupled  mode 
theory,  effective  index,  beam  propagation  method  and  finite-element  analysis,  have  been 
used  [15, 24, 25,  86,  88, 89, 94, 100-125],  Analysis  of  these  methods  shows  good 
agreement  between  all  of  the  different  models,  provided  that  we  are  far  from  cutoff.  This 
means  that  the  normalized  frequency,  defined  as 

B  =  (bk/7t)[n  -  n(l  -  A)]1/2,  (91) 

is  greater  than  or  equal  to  1 .6  for  single  mode  operation  [114,  125],  In  Eq.  (91),  the  free- 
space  propagation  constant  of  plane  waves  is  defined  as  k  =  2k/a,  or  wave  number,  where  A. 
is  the  wavelength.  Since  the  waveguides  fabricated  for  this  experiment  will  be  rectangular 
ridge  type  structures,  we  will  use  several  of  the  dielectric  strip  waveguide  methods  to  model 
the  devices  for  this  experiment  [25,  88,  89,  102-109]. 

5.3.4  Rectangular  Dielectric  Waveguide  Model 

In  1969,  Marcatili  [88]  presented  analytical  models  for  integrated  optical  rectangular 
dielectric  waveguides  and  directional  couplers.  Fig.  15  illustrates  a  directional  coupler  with 
the  following  parameters:  nj  and  n2  are  the  refractive  indices  of  the  core  and  top  clad, 
respectively;  n3  is  the  refractive  index  of  the  outside  clad;  n4  is  the  refractive  index  of  the 
bottom  clad;  n5  is  the  refractive  index  of  the  center  clad;  the  width  and  thickness  of  the 
individual  waveguides  are  a  and  b,  respectively;  and  the  distance  of  separation  between  the 
two  waveguides  is  c.  The  modes  of  the  dielectric  waveguide  are  classified  as  transverse 
electric  (TEpq)  and  transverse  magnetic  (TMpq)  modes.  The  variables,  p  =  (0,1,2,...)  and  q  = 
(0,1,2,...),  indicate  the  number  of  extrema  of  the  electric  or  magnetic  field  along  the  x  and  z 
directions,  respectively.  This  gives  the  order  of  the  mode.  For  the  fundamental  "0th  order" 
TEpq  mode,  p  =  q  =  0,  and  the  transverse  electric  mode  is  given  the  designation  TEoo.  TEpq 
(TMpq)  modes  have  the  electric  field  component  polarized  predominantly  in  the  z  (x) 
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direction  and  the  magnetic  field  component  polarized  predominantly  in  the  x  (z)  direction. 
The  amplitudes  of  the  fields  have  only  one  maximum  in  each  direction. 


Z 


y6— ►  x 


Figure  15.  Rectangular  Waveguide  Directional  Coupler  Cross  Section. 


For  well-guided  modes,  the  field  decays  exponentially  in  the  n2,  n3,  n4  and  n5  regions 
of  the  directional  coupler  waveguide  structure.  Most  of  the  power  travels  in  the  ni  regions 
while  only  a  small  part  of  the  power  travels  in  the  n2,  n3,  n4  and  n5  regions.  An  even  smaller 
portion  of  the  power  travels  in  the  shaded  areas.  Marcatili  found  that  very  little  error  was 
introduced  when  the  fields  along  the  edges  of  the  shaded  regions  were  not  properly  matched 
when  calculating  the  fields  in  the  ni  regions.  This  allowed  him  to  neglect  the  shaded  areas 
from  consideration,  thus,  simplifying  the  boundary  value  problem  by  solving  only  along  the 
four  sides  of  the  core  region  n^ 
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5.3.5  Marcatili  TEpq  Modes 

For  the  rectangular  dielectric  waveguide,  illustrated  in  Fig.  15,  the  electric  and 
magnetic  field  components  in  the  five  areas  are  given  by  [88] 


Hxi  =  Micos(kxx  +  a)cos(kzz  +  8)exp(-ikyy  +  icot). 

(92) 

Hx2  =  M2cos(kxx  +  a)exp[-(ikz2Z  +  ikyy  -  icot)], 

(93) 

Hx3  =  M3Cos(kzz  +  5)exp[-(ikx3x  +  ikyy  -  icot)], 

(94) 

Hx4  =  M4cos(kxx  +  a)exp(ikz4z  -  ikyy  +  icot), 

(95) 

HX5  =  M5cos(kzz  +  S)sin(kx5x  +  y)exp(-ikyy  +  icot), 

(96) 

O 

1! 

II 

II 

II 

S3 

II 

N 

(97) 

Hyl  =  -[(i/ky)9  2Hx]/9  x5  z]. 

(98) 

Hy2  =  -[(i/ky)  a2Hx2/ax3z], 

(99) 

Hy3  =  -[(i/ky)  d  2Hx3/3x3z], 

(100) 

Hy4  =  -[(i^cy)  a2Hx4/8xaz], 

(101) 

Hy5  =  -[(i/ky)  a  2Hx5/  a  x  a  z] , 

(102) 

Exi  =  -[(l/(008ni2ky))  a  2Hxl/a  xa  z], 

(103) 

Ex2  =  -[(l/(ojen22ky))  a ^/dxdz]. 

(104) 

Ex3  =  -[(l/(C08n32ky))  3  2Hx3/a  x3z], 

(105) 

Ex4  =  -[(l/(C0£n42ky))  a  2Hx4/a  x3z]. 

(106) 

Ex5  =  -[( l/(coen52ky))  a  2Hx5/a  x  d  z]. 

(107) 

Ezi  =  Hxi[(k2n!2  -  kz]2)/coenj2ky)], 

(108) 

Ez2  =  HX2[(k2n22  -  kz22)/co£n22ky)], 

(109) 

Ez3  =  Hx3[(k2n32  -  kz32)/co6n32ky)], 

(HO) 

Ez4  =  Hx4[(k“n4  -  kz42)/co£n42ky)], 

(HI) 

Ezs  =  Hx5[(k2n52  -  k252)/coen52ky)], 

(112) 

Eyi  =  [(i/coeni2)  d  Hxl/a  z]. 

(113) 

Ey2  =  [(i/coen22)  a  Hx2/a  z]. 

(114) 

Ey3  =  [(i/co£n32)aHx3/az], 

(115) 

Ey4  =  [(i/cosn42)  a  Hx4/a  z], 

(116) 
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Ey5  =  [(i/coen52)  d  Hx5/9  z], 


(117) 


where  Mi, 23,4,5  are  the  amplitudes  of  the  field  in  each  respective  medium,  a  locates  the  field 
maxima  in  the  core  region  nj,  8  locates  the  field  minima  in  the  core  region  ni,  co  is  the 
angular  frequency  and  e  is  the  permittivity  of  free  space.  The  y  term  is  necessary  because 
single  mode  operation  of  a  directional  coupler  will  guide  two  kinds  of  TEoo  modes.  One  is 
symmetric  when  y=  90°,  and  the  other  is  antisymmetric  when  y = 0°  [88].  Both  are  TEM 
modes  with  main  field  components  Ey  and  Hx.  For  a  symmetric  mode,  the  plane  x  =  0  is  a 
magnetic  short  circuit.  For  an  antisymmetric  mode,  the  plane  x  =  0  is  an  electric  short 
circuit.  The  electric  and  magnetic  field  components  for  both  symmetric  and  antisymmetric 
modes  are  depicted  in  Fig.  16. 

The  propagation  constants  kxv,  ky  and  kzv  in  each  of  the  five  regions  are  related  by 

[77] 

kxv2  +  ky2  +  kzv2  =  (D2£(xnv2  =  kv2,  (118) 

where  v  =  1, 2, 3, 4, 5  indicates  the  region  and  |_i  is  the  permeability  of  free  space.  In  order 
to  match  the  fields  along  the  boundaries  between  the  core  region  ni  and  the  top  and  bottom 
clad  regions  (m  and  at)  and  between  the  core  region  ni  and  the  outside  and  center  clad 
regions  (n3  and  ns),  Marcatili  assumed  the  following: 


kxi  —  kx?  —  kX4  —  kx 

(119) 

ti 

js 

ii 

11 

3 

(120) 
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Figure  16.  Field  Intensity  Profiles  of  the  Symmetric  and  Antisymmetric  Modes  in  a 
Rectangular  Directional  Coupler. 

Since  the  refractive  index  of  the  core  region  nj  is  larger  than  the  refractive  indices  of 
the  cladding  regions,  only  modes  of  plane  wavelets  impinging  at  grazing  angles  on  the 
surface  of  the  core  region  ni  are  guided.  This  implies  that  kx  «  ky  and  kz «  ky  and, 
therefore,  the  electric  field  components  Exl,  Ex2,  Ex3,  Ex4  and  Ex5  can  be  neglected. 
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Matching  the  remaining  field  components,  Hx  and  Ez,  along  the  boundaries  of  the  core 
region  ni  gives  the  transcendental  equations  [88] 

kxa  =  kx0a[  1  +  (2^/a)exp(-c/^  -  i2y)/(  1  +  kx02^2)]  (121) 

and 

kzb  =  (q  +  l)7t  -  tan‘1[(n22/ni2)(kzri2)3  -  tan'1[(n42/ni2)(kzri4)],  (122) 

where  the  transverse  propagation  constants  kx  and  kz  are  solutions  to  Eqs.  (121)  and  (122), 
respectively.  kxo  is  the  solution  to  [88] 

kx0a  =  (p  +  l)7t  -  tan'*(kxo^3)  -  tan''(kx0^5).  (123) 

A  measure  of  the  penetration  depths  of  the  field  components  in  the  various  media  are  given 
by  [88] 


^3,5=l/[(7t/A3,5)2-kxo2]1/2 

(124) 

and 

Tl2,4=l/[(7t/A2,4)2-kz2]1/2, 

(125) 

where 

A2,3,4,5  =  X/[2(n]2  -  n23A52)1/2] 

(126) 

are  the  maximum  thicknesses  of  the  media  for  which  the  slab  supports  only  the  fundamental 

TEqo  mode. 

From  Eqs.  (1 18)  -  (120),  the  axial  propagation  constant  for  each  mode  of  the  coupler 

is  determined  by  [88] 

ky=(k,2-kx2-kz2)1/2, 

(127) 
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where  ki  =  kni  is  the  propagation  constant  of  a  plane  wave  through  a  medium  with  a 
refractive  index  nj  at  a  free-space  wavelength  X,  and  kx  and  kz  are  obtained  from  Eqs.  (121) 
and  (122). 

Substituting  y  into  Eq.  (121)  gives  slightly  different  axial  propagation  constants  for 
the  symmetric  (y  =  90°)  modes  and  antisymmetric  (y  =  0°)  modes  [88], 

kys  =  kyo[l  +  2(kx02/ky02)(^5/a)exp(-c/^5)/(l  +  kx02^2)] 
and 

kya  =  ky0[  1  -  2(kx02/ky02)(^a)exp(-c/^5)/(l  +  kx02^2)], 

where 

ky0  =  (kj2  -  kx02  -  kz2)1/2  (130) 

is  the  axial  propagation  constant  of  the  TE  mode  of  a  single  dielectric  waveguide  (c  — >  ) 

(see  Fig.  16.). 

The  interaction  length  L  for  total  power  transfer  from  one  waveguide  to  the  other 
and  the  coupling  coefficient  K  between  the  two  waveguides  can  then  be  related  to  kys  and 
kya  by  [88] 

-iK  =  ti/(2L)  =  (kyS  -  kya)/2  =  2(kx02/kyo)(^a)exp(-c/^)/(l  +  kx02^2).  (131) 

The  interaction  length  and  coupling  coefficient  for  the  zero-gap  directional  coupler  are 
obtained  by  driving  the  separation  between  the  two  adjacent  waveguides  to  zero  and 
substituting  the  refractive  index  n4  for  ns,  which  yields  [32] 

-iK  =  n/(2L)  =  2(kx02/kyo)£s/a)/(l  +  kx0V).  (132) 


(128) 

(129) 
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5.3.6  Marcatili  TMpq  Modes 

The  analysis  for  the  TMpq  modes  is  very  similar  to  that  used  for  the  TEpq  modes. 

The  electric  and  magnetic  field  components  for  TMpq  modes  are  found  by  substituting  E  for 
H,  H  for  E,  -£  for  ji  and  -ji  for  £  into  Eqs.  (92)  -  (1 17).  Following  Eqs.  (119)  and  (120),  the 
field  components  E&  and  Hxv  can  be  neglected  for  the  TMpq  analysis. 

Matching  the  field  components,  Ex  and  Hz,  along  the  boundaries  of  the  core  region 
ni  in  the  same  way  as  was  done  to  match  the  field  components  Hx  and  Ez  for  the  TEpq 
analysis  gives  the  transcendental  equations  [88] 

kxa  =  kx0a[  1  +  (2^/a)exp(-c/^  -  2iy)/(  1  +  kx02^2)]  (133) 

and 

kzb  =  (q  +  l)7t  -  tan''(kzri2)  -  tan_1(kzr|4),  (134) 

where  the  transverse  propagation  constants  kx  and  kz  are  solutions  to  Eqs.  (133)  and  (134), 
respectively.  kxo  is  the  solution  to 

kx0a  =  (p  +  1)7T  -  tan'1[(n32/ni2)(kx0^3)]  -  tan'1[(n52/n12)(kxo^5)].  (135) 

r|2, 114,  ^3  and  ^  are  obtained  from  Eqs.  (124)  -  (126). 

The  interaction  length  L  for  total  power  transfer  from  one  waveguide  to  the  other 
and  the  coupling  coefficient  K  between  the  two  waveguides  are  determined  by  Eq.  (131) 
using  kyo,  kz  and  kx0  obtained  from  Eqs.  (130),  (134)  and  (135),  respectively.  The 
interaction  length  and  coupling  coefficient  for  the  zero-gap  directional  coupler  are  obtained 
by  substituting  ky0,  kz  and  kxo,  obtained  from  Eqs.  (130),  (134)  and  (135),  into  Eq.  (132) 

[32]. 


Marcuse  [89]  derived  a  rectangular  dielectric  waveguide  model  following  the 
approach  developed  by  Marcatili.  It  is  also  widely  used  and  gives  identical  results. 


55 


5.3.7  Effective  Index  Method 


In  a  1975  paper,  Bums  and  Milton  discussed  mode  conversion  in  planar  dielectric 
separating  waveguides  [102],  This  has  since  become  known  as  the  effective  index  method, 
and  is  treated  in  detail  by  Nishihara,  Haruna  and  Suhara  for  both  two  and  three-  dimensional 
dielectric  waveguide  structures  [25],  The  effective  index  method  is  based  on  the  concept  of 
the  effective  index  of  the  guided  mode.  The  guided  mode  propagating  along  the  z  axis  sees 
the  effective  index  (N)  defined  by  [25, 121] 

P  =  kN>  (136) 

where  (3  is  the  axial  propagation  constant,  k  is  the  free  space  propagation  constant,  and  N  = 
nicos0i  (see  Fig.  17). 


Guided  Wave 


z 

▲ 

y  6 - ►  X 

Figure  17.  Trajectory  of  a  Guided  Wave  in  a  Slab  Waveguide. 
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The  effective  index  method  used  for  the  dielectric  waveguide  structure  fabricated  for 
this  experiment  (illustrated  in  Fig.  15)  is  performed  by  dividing  the  three-dimensional 
waveguide  into  two  two-dimensional  waveguides.  The  first  is  2-D  Waveguide  1,  with  light 
confinement  in  the  z  direction.  The  second  is  2-D  Waveguide  2,  with  light  confinement  in 
the  x  direction  (see  Fig.  18).  The  axial  propagation  constant  for  the  TEN  modes  of  a  three- 
dimensional  waveguide  structure  is  determined  using  the  following  procedure:  The 
effective  index  Ni  of  2-D  Waveguide  1  is  found  by  applying  TMpq  mode  analysis  to  the  top 
and  bottom  clad  regions  (n2  and  n 4,  respectively)  with  waveguide  thickness  b.  Effective 
index  Ni  is  used  to  find  the  effective  index  N2  of  2-D  Waveguide  2  by  applying  TEpq  mode 
analysis  to  the  outside  and  center  clad  regions  (n3  and  n5,  respectively)  with  waveguide 
width  a.  The  axial  propagation  constant  is  obtained  by  substituting  N2  for  N  in  Eq.  (136). 


Three-Dimensional  Dielectric  Waveguide 


n4 


Two-Dimensional  Waveguide  1 


Two-Dimensional  Waveguide  2 


Figure  18.  Analytical  Model  for  the  Effective  Index  Method  for  TEpq  Modes. 
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The  axial  propagation  constant  for  the  TMpq  modes  is  determined  in  a  manner 
similar  to  the  one  used  for  the  TEpq  modes.  Effective  index  Nt  of  2-D  Waveguide  1  is  found 
by  applying  TE^  mode  analysis  to  the  top  and  bottom  clad  regions  (n2  and  m,  respectively) 
with  waveguide  thickness  b.  Effective  index  Ni  is  then  used  to  find  effective  index  N2  of  a 
2-D  Waveguide  2  by  applying  TMpq  mode  analysis  to  the  outside  and  center  clad  regions  (n3 
and  n5,  respectively)  with  waveguide  width  a. 

5.3.8  Effective  Index  Method  TEpq  Modes 


As  outlined  in  the  preceding  section,  2-D  Waveguide  1  represents  the  TMpq  modes 
in  an  infinite  slab  waveguide  of  thickness  b.  2-D  Waveguide  2  represents  the  TEpq  modes  in 
an  infinite  slab  waveguide  of  width  a.  Since  Ey  is  the  dominant  electric  field  component,  the 
field  appears  as  a  TEpq  mode  with  respect  to  the  outside  and  center  clad  regions  (n3  and  n5, 
respectively).  The  analysis  can  be  represented  by  the  following  [1 14]: 

P2  =  k[N2{TMb(n2ln1ln4);TEa(n3IN,ln5)}],  (137) 

where  (32  is  the  TEpq  axial  propagation  constant  for  the  three  dimensional  dielectric 
waveguide. 

First,  we  start  with  the  wave  equation  for  the  TMpq  modes  of  2-D  Waveguide  1 , 
given  by  [25] 


a2Hx/3z2  +  (k2n2-p,2)Hx  =  0, 


(138) 


where 

Ez  =  (pi/G0£on2)Hx, 
and 

Ey  =  ( l/ico£0n2)  d  Hx/ 5  z, 

where  is  the  TMpq  axial  propagation  constant  for  2-D  waveguide  1. 


(139) 

(140) 
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The  field  solutions  can  then  be  written  as  [25] 


Hxl  =  HiCos(kzz  +  (jb),  (141) 

at  -b  <  z  <  0  for  the  core  region  ni, 

=  H2exp(-y2z),  (142) 

at  z  >  0  for  the  top  clad  region  n2,  and 

Hz4  =  H4exp[y4(z  +  b)] ,  (1 43) 

at  z  <  -b  for  the  bottom  clad  region  n4.  The  variables  y2  and  y4  are  the  transverse 
propagation  constants  of  a  plane  wave  in  the  top  and  bottom  cladding  regions,  respectively. 
<t>2  is  the  phase  shift  suffered  by  a  wave  polarized  with  its  magnetic  field  vector  parallel  to 
the  interface  between  the  core  region  nj  and  the  top  clad  region  n2. 

In  addition  to  the  phase  shift  due  to  reflection  of  the  wave  from  the  interface 
between  ni  and  n2,  there  is  also  a  phase  shift  <|)4  suffered  by  the  reflection  of  the  wave  from 
the  interface  between  the  core  region  ni  and  the  bottom  clad  region  n4.  Using  Eqs.  (138)  - 
(143) ,  these  phase  shifts  can  be  defined  mathematically  by  [89] 

((h  =  -2arctan[(ni2/n22)(p!2  -  n22k2)1/2/(ni2k2  -  Pi2)1/2]  (144) 

and 

(j)4  =  -2arctan[(ni2/n42)(Pi2-n42k2)I/2/(ni2k2- pj2)1/2].  (145) 

Using  the  differences  between  the  phase  fronts,  which  are  due  to  the  reflection  of  the 
wave  at  the  two  interfaces  and  given  mathematically  by  [89] 

Si  =  (bcos0i)[(l/tan0i)  -  tan0i]  =  (b/sin0O(cos20i  -  sin20O  (146) 

and 

s2  =  b/sin0j,  (147) 
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the  condition  in  which  both  reflected  waves  contribute  to  the  same  plane  waves  can  be 
expressed  by  the  relation  [89] 

ni(s2  -  sj)k  +  <t>2  +  <j)4  =  2(q  +  \)k,  (148) 

where  q  =  (0,1,2,...)  indicates  the  number  of  extrema  of  the  magnetic  field  in  the  y  direction. 

The  eigenvalue  equation,  obtained  from  Eqs.  (144)  -  (148),  can  then  be  written  as 
[25] 


kzb  =  (q  +  l)7t  -  tan‘1[(n22/n12)(k2/Y2)]  -  tan'1[(n42/n12)(kz/Y4)],  (149) 


where 

kz  =  (n,2k2-p,2)1/2, 


(150) 


2i,2\ 1/2 


Y2  =  ((V -n2V) 


and 


Y4  =  (Pi2  -  n42k2)1/2 


(151) 

(152) 


are  the  transverse  propagation  constants. 

Using  Eq.  (136),  the  transverse  propagation  constants,  given  by  Eqs.  (150)  -  (152), 
can  be  expressed  in  terms  of  the  effective  index  of  2-D  Waveguide  1 . 


kz  =  k(n,2-N12)1/2, 

(153) 

Y2  =  k(Ni2-n22)1/2 

(154) 

ii 

2 

1 

I5 

^JO 

to 

(155) 
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Next,  the  wave  equation  for  the  TEpq  modes  of  2-D  Waveguide  2  is  given  by  [25] 
92Ez/ax2  +  (k2n2-p22)Ez  =  0,  (156) 


where 


Hx  =  -(p2/cofio)Ez 

(157) 

and 

Hy  =  -(1/icojj.o)  d  Ez/  d  x. 

(158) 

The  field  solutions  can  be  written  as  [25] 

Ezl  =  E!Cos(kxx  +  <t>3), 

(159) 

at  0  <  x  <  a  for  the  core  region  now  defined  as  Nj, 

Erf  =  E2exp(-y3x), 

(160) 

at  x  <  0  for  the  outside  clad  region  n3,  and 

Ezs  =  E5exp[7s(x  +  a)]. 

(161) 

at  x  >  a  for  the  center  clad  region  ns.  The  variables  73  and  75  are  the  transverse  propagation 
constants  of  a  plane  wave  in  the  outside  and  center  cladding  regions,  respectively,  (fe  is  the 
phase  shift  suffered  by  a  wave  polarized  with  its  electric  field  vector  parallel  to  the  interface 
between  the  core  region  Nj  and  the  outside  clad  region  n3. 

In  addition  to  the  phase  shift  due  to  reflection  of  the  wave  from  the  interface 
between  Ni  and  n3,  there  is  also  a  phase  shift  tps  suffered  by  the  reflection  of  the  wave  from 
the  interface  between  the  core  region  Nj  and  the  center  clad  region  ns.  Using  Eqs.  (156)  - 


61 


(161),  these  phase  shifts  can  be  defined  mathematically  by  [25] 


<b  =  -2arctan[(P22  -  n32k2)1/2/(ni2k2  -  p22)1/2]  (162)  . 

and 

(J)5  =  -2arctan[(p22  -  n52k2)1/2/(ni2k2  -  p22)172].  (163) 

Using  the  differences  between  the  phase  fronts,  which  are  due  to  the  reflection  of  the 
wave  at  the  two  interfaces  and  given  mathematically  by  [25] 

Si  =  (acos0i)[(l/tan9i)  -  tan0i]  =  (a/sin0i)(cos20i  -  sin20O  (164) 

and 

s2  =  a/sin0,,  (165) 

the  condition  in  which  both  reflected  waves  contribute  to  the  same  plane  waves  can  be 
expressed  with  the  relation  [25] 


ni(s2  -  Si)k  +  4>3  +  (j>5  =  2(p  +  1)tt,  (166) 

where  p  =  (0,1,2,...)  indicates  the  number  of  extrema  of  the  electric  field  in  the  x  direction. 

The  eigenvalue  equation,  obtained  from  Eqs.  (162)  -  (166),  can  then  be  written  as 
[25] 


kxa  =  (p  +  l)7t  -  tan‘'(kx/y3)  -  tan'Vk/ys),  (167) 

where 


kx  =  (N,2k2-|322)1/2, 

(168) 

Y3  =  (P22-n32k2)I/2 

(169) 

y5  =  (p22-n52k2),/2 

(170) 
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are  the  transverse  propagation  constants. 


Using  Eq.  (136),  the  transverse  propagation  constants,  given  by  Eqs.  (168)  -  (170), 
can  be  expressed  in  terms  of  the  effective  index  of  2-D  Waveguide  2. 

kz  =  k(N,2-N22)1/2,  (171) 

Y3  =  k(N22-n32)1/2  (172) 

and 

Y5  =  (N22-n52)1/2.  (173) 

Having  solved  for  p2  and  kx,  the  symmetric  propagation  constant  ps  and 
antisymmetric  propagation  constant  pa  (corresponding  to  the  symmetric  and  antisymmetric 
TEoo  modes  of  the  directional  coupler,  respectively)  are  found  by  [88] 

Ps  =  P2[l  +  2(kx2/p22)(l/aY5)exp(-cY5)/(l  +  kx  V)] 
and 

Pa  =  p2[l  -  2(kx2/p22)(l/aY5)exp(-cY5)/(l  +  kxV)L 
where  c  is  the  distance  of  separation  between  the  two  waveguides. 

The  interaction  length  L  for  total  power  transfer  from  one  waveguide  to  the  other, 
and  the  coupling  coefficient  K  between  the  two  waveguides,  can  be  related  to  P$  and  pa  by 
[25] 

-iK  =  rt/(2L)  =  (Ps  -  Pa)/2.  (176) 

The  interaction  length  and  coupling  coefficient  for  the  zero-gap  directional  coupler 
are  obtained  by  driving  the  separation  between  the  two  adjacent  waveguides  to  zero  and 


(174) 

(175) 
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substituting  the  refractive  index  n4  for  n5  [32].  This  yields, 

(177) 

(178) 

which  are  then  substituted  back  into  Eq.  (176). 


(3S  =  p2[l  +  2(kx2/p22)(l/ay5)(l/(l  +  kx2/y52))] 
and 

|3a  =  p2[l  -  2(kx2/p22)(l/ay5)(l/(l  +  kx2/y52))]. 


5.3.9  Effective  Index  Method  TMN  Modes 


For  TMpq  mode  analysis,  2-D  Waveguide  1  represents  the  TEpq  modes  in  an  infinite 
slab  waveguide  of  thickness  b.  2-D  Waveguide  2  represents  the  TMpq  modes  in  an  infinite 
slab  waveguide  of  width  a.  With  Hy  as  the  dominant  magnetic  field  component,  the  field 
appears  as  a  TMpq  mode  with  respect  to  the  outside  and  center  clad  regions  (n3  and  n5).  The 
analysis  can  be  represented  by  the  following  [114]: 

P2  =  k[N2{TEb(n2lniln4);TMa(n3lN,ln5)}],  (179) 

where  p2  is  the  TMpq  axial  propagation  constant  for  the  three  dimensional  dielectric 
waveguide.  Following  the  same  procedure  as  was  presented  in  the  previous  section,  the 
eigenvalue  equation  for  the  TEpq  modes  of  2-D  Waveguide  1  can  be  written  as  [25] 

kxa  =  (p  +  l)7t  -  tan‘1[(n32/ni2)(kx/y3)]  -  tan''[(n52/ni2)(kx/y5)],  (180) 


where 

kx  =  (nfk2  -  p,2)1/2, 


(181) 


and 


73  =  (Pi2  -  n32k2)1/2 
y5  =  (p,2-n52k2)1/2. 


(182) 

(183) 
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Using  Eq.  (136),  the  transverse  propagation  constants,  given  by  Eqs.  (181)  -  (183), 
can  be  expressed  in  terms  of  the  effective  index  of  2-D  Waveguide  1. 


kz  =  k(n12-N12)1/2, 


(184) 


Y3  =  k(N,2-n32)1/2 


(185) 


y5  =  (N12-n52)1/2. 


(186) 


The  eigenvalue  equation  for  the  TMpq  modes  of  2-D  Waveguide  2  can  be  written  as  [25] 


kzb  =  (q  +  l)7t  -  tan''(kz/y2)  -  tan'’(kz/y4), 


(187) 


where 


kz  =  (N,2k2-p22)1/2 


(188) 


y2  =  (Pi2-n22k2)1/2 


(189) 


y4  =  (P22  -  afk2)172. 


(190) 


Using  Eq.  (136),  the  transverse  propagation  constants,  given  by  Eqs.  (188)  -  (190), 
can  be  expressed  in  terms  of  the  effective  index  of  2-D  Waveguide  2. 


kz  =  k(N]2-N22)1/2, 


(191) 


y2  =  k(N22-n22)1/2 


(192) 


/XT  2  2\l/2 

y4  =  (N2-  -  n4-)i/". 


(193) 
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Having  solved  for  p2  and  kx,  the  symmetric  propagation  constant  (3S  and  antisymmetric 
propagation  constant  |3a  for  the  TMoo  modes  of  the  directional  coupler  are  found  by 
substituting  fc  and  kx,  determined  here,  into  Eqs.  (174)  and  (175). 

The  interaction  length  L  for  total  power  transfer  from  one  waveguide  to  the  other, 
and  the  coupling  coefficient  K  between  the  two  waveguides,  can  be  related  to  (3S  and  pa  by 
Eq.  (176).  The  interaction  length  and  coupling  coefficient  for  the  zero-gap  directional 
coupler  are  obtained  by  using  Eqs.  (177)  and  (178),  which  are  then  substituted  back  into  the 
equation  for  the  zero-gap  directional  coupler,  Eq.  (176). 

5.3.10  Beam  Propagation  Method 

In  1976,  Fleck,  Morris  and  Feit  proposed  a  propagating  beam  method  (bpm)  for 
computing  fields  of  an  optical  beam  passing  through  a  medium  [114].  The  method 
calculates  the  characteristics  of  an  input  beam  as  it  propagates  through  a  medium  over  a 
small  distance  y,  and  then  corrects  for  the  Variations  of  the  refractive  index  as  seen  by  the 
beam  over  this  distance  y.  It  has  already  been  applied  to  several  optical  fiber  and  integrated 
optic  applications  [104-110, 117,  123, 124], 

The  beam  propagation  method  is  derived  for  a  scalar  field.  The  theory  is  restricted 
to  small  changes  in  the  refractive  index.  The  first  part  of  the  derivation  assumes  the 
propagation  of  a  high  frequency  beam  through  an  inhomogeneous  medium.  It  begins  with 
the  wave  equation  [109] 

V  2(j)  +  k2n2(r)(|)  =  0,  (1 94) 

where  (j>  represents  the  scalar  field,  k  is  the  free  space  wave  number  defined  at  the  beginning 
of  this  chapter,  and  n(r)  =  n0(r)2  +  An2(r)  is  the  refractive  index  of  the  medium  consisting  of 
an  unperturbed  part  (n0(r)2),  and  a  perturbation  (An2(r)).  Splitting  n2(r)  into  its  respective 
parts,  Eq.  (194)  takes  the  form  [109] 

V  24>  +  k2n02(r)(j>  =  -K2Av2(r)(J)  =  s(r),  ( 1 95) 


66 


where  the  right  hand  side  of  the  equation,  designated  as  s(r),  is  considered  a  source  function. 
Next,  the  unperturbed  part  of  the  index  is  chosen  so  that 
the  new  wave  equation  [109] 

V  +  k2no2(r)\|/  =  0,  (1 96) 

along  with  the  radiation  conditions  at  y  =  has  a  solution.  If  V|/  is  known  for  y  =  yo,  where 
yo  is  the  position  y  =  0,  then  \jr  and  its  derivative,  with  respect  to  y,  can  be  solved  for  all 
values  of  y  by  use  of  an  operator  (a),  a  acts  with  respect  to  only  the  transverse  coordinates 
(x,z).  This  leads  to  [109] 

d\\f/dz  =  a\|/(x,y0,z).  (197) 

For  n0  =  constant,  the  scalar  field  can  be  represented  in  the  angular  spectrum  domain  by 
[109, 126] 


'F  (kx,kz,y)  =  JJ  \|/(x,y,z)exp[-i(kxx  +  kzz)]dxdz, 


(198) 


where  TTk^k^y)  is  the  Fourier  transform  of  the  scalar  field  (iy(x,y,z)),  and  kx  and  kz  are  the 
transverse  propagation  constants.  Taking  the  derivative  of  Eq.  (198),  with  respect  to  y, 
yields  [109, 126] 

a  TV3  y(kx,k2,y)  =  -iky^(kx,k2,y),  (199) 

where  ky  =  (k2no2  -  kx2  -  kz2)1/2  is  the  axial  propagation  constant.  By  using  Eq.  (199),  the 
integral  operator  a  can  be  derived. 

For  a  given  coordinate  (y),  the  field  (()>)  is  split  into  two  parts.  The  first  part,  (<j>i),  is 
generated  by  the  sources  in  the  region  where  y'  <  y.  The  second  part,  (<jb),  is  generated  by 
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the  sources  where  y'  >  y.  Using  the  function  [109] 


(ei(rlr')  =  0  for  y  <  y' 
ei(rlr')  =  1/2  for  y  =  y' 

ei(rlr0=  1  fory>y'},  (200) 

an  explicit  expression  for  ^  and  fa  can  be  obtained.  Defining  G  as  Greens  function  of  Eq. 
(196),  the  wave  equation  (fo)  can  be  expressed  by  [109] 


+CO 

<th(r)  =  J  JJ  G(rlrOei(rlrOs(r')dV'. 


Taking  the  derivative  of  (j)h  with  respect  to  y,  yields  [109] 


9  <|>i/ 9  y(r)  =  J  J  J  9  G/9  yCrlrOedrlrOsOOdV' 

— OO 

+oo 

+  J  JJ  G(rlr')5(y  -  y0s(r9dV" 


(201) 


(202) 


where  5(y  -  y  )  is  the  step  size.  The  first  integral  of  Eq.  (202)  represents  the  unperturbed 
medium.  It  can  be  expressed  in  terms  of  a  using  Eq.  (197)  [109], 

a<h(r).  (203) 

The  second  integral  of  Eq.  (202)  includes  both  the  unperturbed  medium  and  the 
perturbation.  With  the  help  of  Eq.  (195),  the  second  integral  in  Eq.  (202)  can  be  written  as 
[109] 


+oo 

JJ  G(x,y,zlx',y,z9(-k2)An2(x',y,z')())(x',y,z9dx'dz".  (204) 
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Eq.  (204)  can  be  expressed  in  terms  of  an  operator  (u)  acting  on  (f)  by  [109] 


u<|)(x',y,zO,  (205) 

where,  like  a,  operator  u  acts  with  respect  to  only  the  transverse  coordinates  (x,z). 
Substitution  of  Eqs.  (203)  and  (205)  into  Eq.  (202)  yields  [109] 

d  <>i/d  y(x",y,z")  =  a<|)i(x',y,z')  +  fi<|)(x',y,z').  (206) 

Neglecting  the  influence  of  the  reflected  field  (t^)  on  <J>i,  4>i  can  be  substituted  for  <j)  in  Eq. 
(206).  The  first-order  differential  equation  for  4>i  then  becomes  [109] 

d  fa/d  y(x',y,z')  =  a<(>i(x',y,z')  +  u<t>i(x',y,zO,  (207) 

where  u<}>i  is  the  correction  term. 

From  Eq.  (207),  the  propagation  of  field  (4>i)  is  given  by  two  terms.  One  describes 
propagation  in  an  unperturbed  medium  (a<t>i).  The  other  describes  a  correction  term  (u<|)i), 
representing  the  influence  of  An.  Using  Fourier  transforms  and  convolution,  along  with  the 
assumption  that  the  unperturbed  media  is  homogeneous,  the  correction  term  in  Eq.  (207)  can 
be  simplified  [109] 

fl<t>i(x',y,zO  =  -(ik/2n0)An2(|>i(x',y,z')-  (208) 

By  combining  Eqs.  (207)  and  (208),  the  first-order  differential  equation  for  <j>j  becomes 
[109] 


d  <j>i/d  y(x',y,zO  =  a<{>,(x',y,z')  -  (ik/2n0)An2<t>i(x',y,z').  (209) 

In  order  to  solve  Eq.  (209),  it  will  be  assumed  that  <J>i  is  known  at  position  y  =  yo.  If  An2  =  0, 
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Eq.  (209)  reduces  to  [98] 


a  e/a  y(x,y0,z)  =  ae(x,y0,z),  (210) 

where  e(x,yo,z)  =  <()i(x,yo,z)  is  given.  Therefore,  e  represents  the  field  as  it  propagates 
through  a  medium  with  a  refractive  index  (n0(r)).  Introducing  a  correction  factor  (exp(T)) 
leads  to  an  approximate  solution  to  the  differential  Eq.  (209).  This  solution  is  derived  by 
[109] 


<t>i(r)  =  exp[r(r)]e(r),  (211) 

with  T(x,y0,z)  =  0.  Combining  Eqs.  (209)  and  (21 1)  yields  [109] 

a  17  3  y  =  -(ik/2n0)An2  +  [a(exp(T)e)  -  exp(T)a(e)]/exp(r)e-  (2 1 2) 

It  has  been  demonstrated  that  the  second  term  in  Eq.  (212)  can  be  neglected,  provided  Eq. 
(212)  is  integrated  over  a  small  distance  (y  -  y0)  [109].  This  yields 

T  =  -(ik/2n0)An2(y  -  y0),  (2 1 3) 

for  (y  -  y0)  small.  Combining  Eqs.  (211)  and  (213)  yields  [109] 

<|)i(x,yo,z  +  Ay)  =  £(x,y0,z  +  Ay)exp[-(ik/2n0)An2Ay),  (214) 

which  is  the  equation  for  the  correction  term. 

For  n0(r)  independent  of  y,  the  eigenmode  of  Eq.  (196)  can  be  determined.  This 
allows  one  to  calculate  the  propagation  of  £  [109]. 

£(x,y0,z  +  Ay)  =  An\|/n(x,z)exp[-(iKvA\|/)]  (215) 
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where 


-h» 

An  =  jj  £(x,y0,z)\|/„(x,z)dxdz. 


(216) 


Using  Eq.  (215)  and  multiplication  with  the  correction  term,  defined  by  Eq.  (214),  allows 
the  propagation  of  light  in  small  steps. 

Application  of  the  beam  propagation  method  to  media  with  small  refractive  index 
variations  allows  one  to  analyze  beam  propagation  through  integrated  waveguide  structures. 
For  this  case,  the  operator  a  can  be  expressed  in  terms  of  Fourier  transforms  (see  Eqs.  (198) 
and  (199)).  The  propagation  step  defined  by  Eq.  (215)  can,  thus,  be  calculated  using  fast 
Fourier  transforms  (FFT's).  Performing  a  two  dimensional  EFT  extends  the  beam 
propagation  method  to  three  dimensions  (3D),  allowing  one  to  analyze  the  dielectric 
waveguide  structure  illustrated  in  Fig.  15. 

The  beam  propagation  method  for  integrated  optic  waveguide  applications  is 
available  as  computer  aided  design  CAD  software  tools.  It  is  readily  available  for  personal 
computers  (PCs),  Macintosh  (MAC)  computers  and  workstations.  The  BPM  CAD  tool 
utilized  for  this  chapter  was  BeamPROP  for  PC  [127].  Although  the  BPM  does  not 
generate  hard  numbers  for  the  interaction  length  of  directional  couplers,  in  the  way  the 
dielectric  waveguide  and  effective  index  methods  do,  the  beam  propagation  method  does 
provide  an  excellent  prediction  of  device  performance.  It  calculates  the  amplitudes  of  the 
fields  as  they  propagate  through  the  medium.  For  example,  the  BPM  CAD  program  was 
used  to  generate  Fig.  19,  which  gives  a  visual  prediction  of  the  behavior  of  a  light  wave 
propagating  through  an  integrated  optic  waveguide  with  an  S-bend.  The  following 
parameters  were  substituted  into  the  BPM  software:  the  waveguide  width  a  =  3  pm,  the 
waveguide  core  thickness  b  =  1  pm,  the  operating  wavelength  X  =  830  nm,  the  refractive 
index  of  the  core  ni  =  1 .63,  the  refractive  index  of  the  top  and  bottom  clad  regions  n2  =  114  = 
1 .61 ,  and  the  refractive  index  of  the  outside  clad  regions  m  =  ns  =  1.61.  Zero  propagation 
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loss  was  assumed  for  the  straight  guides.  The  S-bend  was  designed  using  the  formula  [128] 


x(y)  =  (h/l)y  -  (h/2ft)sin[(2jt/l)y],  (217) 

where  1  is  the  longitudinal  separation  in  the  y  direction  and  h  is  the  lateral  offset  in  the  x 
direction.  Eq.  (217)  allows  one  to  minimize  the  radiative  loss  and  maximize  the  mode 
confinement.  For  the  S-bend  structure  illustrated  in  Fig.  13,  h  =  20  (im  and  1  =  500  pm. 


Figure  19.  Amplitudes  of  the  Propagating  Fields  in  an  Integrated  Optic  S-Bend  Waveguide. 
5.3.11  Refractive  Index  of  the  Waveguide  Cladding 

The  refractive  index  of  the  cladding  layer  can  be  found  using  the  relation  [16] 
iW  -  ndad  >  m2X2/[4(tg)2(ncore  +  ncIad)],  (2 1 8) 
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where  m  is  the  mode  number  and  tg  is  the  thickness  of  the  core.  Assuming  single  mode 
operation  (m=l),  Eq.  (218)  becomes 


nciad  >  [(ncore)2  -  X2/[4(tg)2]1/2.  (219) 

One  can  also  use  a  TE/TM  mode  analysis  presented  by  Ghatak,  Thyagarajan  and  Shenoy 
[129]  to  model  the  loss  for  a  given  mode. 

The  refractive  indices  and  EO  coefficients  for  the  various  polymers  are  usually 
provided  by  the  manufacturer,  but  they  can  also  be  determined  experimentally.  The 
parameters  used  for  examples  in  this  chapter  are  typical  values.  Now  that  the  basic  theory 
for  designing  EO  devices  has  been  presented,  one  can  address  fabrication  issues. 
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5.4  DEVICE  FABRICATION 


NLO  polymers  typically  come  in  crystal  or  powdered  form.  For  our  example,  we 
will  use  an  NLO  polymer  from  AdTech  Systems  Research  known  as  LD-3  [7],  It  is 
composed  of  a  methymethacrylate  (MMA)  polymer  containing  4-[bis  (2-hydroxyethyl) 
ammo]-4'-  ((6-methacryloyl-hexyl)  sulfonyl)  azobenzene  nonlinear  optical  chromophores.  A 
recipe  for  reconstituting  the  polymer  into  a  liquid  for  spin  casting  is  usually  provided  by  the 
polymer  manufacturer.  We  will  start  with  LD-3  in  powder  form.  First,  grind  0.10  grams  of 
LD-3  into  a  fine  powder  and  dissolve  the  powder  in  1.0- 1.5  ml  of  pure  and  dry  organic 
solvent  tetrahydrofurane  (THF)  [7],  Add  0.060  grams  of  cross-linker  4,4'-diisocyanato- 
3,3'dimethoxybiphenyl,  which  forms  3D  networks  and  locks  the  dipole  aligned 
chromophores  [7].  Stir  the  mixture  for  15  minutes  at  room  temperature.  Filter  the  mixture 
through  a  0.2  pm  syringe  filter  [7].  With  the  material  in  its  liquid  form  we  are  ready  to 
begin  fabrication  of  our  EO  device.  There  are  several  methods  that  have  been  employed  for 
fabricating  NLO  polymer  devices.  The  two  most  common  techniques,  namely  photobleach 
[4, 8,  21, 22,  35,  36, 40, 42, 43, 45, 48,  51-62]  and  etch-and-fill  [8,  21,  22,  35,  39, 45-47, 49- 
54]  are  described  below  in  5.4.1. 

5.4.1  Photobleach 

To  begin  device  processing,  we  evaporate  a  thin  metal  layer  on  top  of  a  substrate, 
which  can  either  be  glass,  polymer  or  semiconductor.  This  metal  layer  will  be  the  ground 
plane  electrode.  It  can  be  gold,  aluminum,  titanium/tungsten/gold  or  any  good  conductor;  it 
should  be  4000  A  -  1  urn  thick  for  good  conductivity.  Deposit  a  passive  polymer  bottom 
cladding  material  with  a  refractive  index  about  0.02-0.04  less  than  that  of  the  core.  The 
cladding  material  is  usually  synthesized  to  work  with  the  NLO  polymer  being  used.  It  can 
be  an  acrylate,  polyimide,  or  polyetherimide,  just  to  name  a  few.  It  must  also  have  a  higher 
glass  transition  temperature  (Tg)  than  the  core  material  in  order  to  perform  the  poling  step. 

It  is  spin  cast  over  the  metal  layer  at  an  acceleration  and  a  speed  to  yield  a  thickness  of  3 
pm.  The  cladding  layer  is  then  cured  by  baking  on  a  hot  plate  until  dry.  The  exact  time 
depends  on  the  type  of  polymer  used.  The  NLO  polymer  core  layer  is  then  deposited  on  top 
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of  the  cladding  layer  at  an  acceleration  and  speed  in  order  to  yield  a  thickness  of  2  |im.  It  is 
baked  until  dry.  A  passive  polymer  top  cladding  layer  is  deposited  in  the  same  manner  in 
order  to  yield  the  same  thickness  as  the  bottom  cladding  layer.  The  top  cladding  layer  is 
also  baked  until  it  is  dry.  A  second  metal  layer,  used  for  waveguide  patterning,  is 
evaporated  on  top  of  the  top  cladding  layer.  It  is  a  good  conductor,  typically  500  A  thick. 
The  metal  layer  is  then  patterned  with  the  desired  EO  waveguide  circuit,  using  standard 
lithographic  techniques,  leaving  a  metal  mask.  The  wafer  is  then  flood  exposed  with 
ultraviolet  (UV)  light.  Exposure  to  UV  results  in  a  slight  and  very  controllable  reduction  in 
the  refractive  index  of  the  NLO  polymer.  This  process  is  called  photobleach. 
Photobleaching  defines  the  side  cladding  layers  of  the  channel  waveguide.  Exposures  take 
several  hours,  and  the  refractive  index  of  the  side  cladding  can  be  controlled  to  the  4th 
decimal  place.  The  metal  mask  used  for  photobleaching  is  then  stripped  off.  A  third  500- 
3000  A  metal  layer  for  the  electrodes  is  evaporated  on  top  of  the  top  clad.  The  top  metal 
layer  is  patterned  with  the  desired  electrode  structure  using  standard  IC  tools  and  techniques. 
After  the  electrodes  have  been  patterned,  the  NLO  material  directly  under  the  electrodes 
needs  to  be  poled  in  order  to  operate  as  an  EO  device.  Poling  is  a  technique  used  to  align 
the  molecules  or  moieties  of  the  material  in  order  to  maximize  the  EO  properties  of  the  NLO 
polymer.  Poling  is  performed  by  heating  the  EO  circuit  to  a  temperature  above  the  glass 
transition  temperature  of  the  material,  specified  by  the  manufacturer,  and  then  applying  a 
DC  voltage.  The  higher  temperature  allows  the  molecules  within  the  material  to  move 
freely.  For  the  NLO  polymer  LD-3,  a  poling  temperature  of  150°C  is  specified.  An  electric 
field  is  then  applied  across  the  electrodes  using  soft  tip  probes.  The  moieties  will  align  in 
the  direction  of  the  electric  field,  thus,  increasing  the  electro-optic  coefficient  r 33  in  those 
regions.  The  poling  voltage  is  typically  150  VDC/pm.  Thus,  for  this  example,  the  poling 
voltage  required  would  be  1200  VDC  [3, 4, 7,  8, 21, 22, 35, 36, 39-42, 48-63, 65].  With  the 
voltage  still  applied,  the  circuit  is  cooled  to  room  temperature.  This  locks  the  moieties  in 
place.  This  type  of  poling  is  known  as  selective  poling  [3, 4,  8,  21, 22,  36, 40, 42, 43, 48, 
51-62];  only  the  NLO  material  under  the  electrode  is  poled.  Since  the  regions  along  the  axis 
parallel  to  the  applied  voltage  become  birefringent,  the  electric  field  of  a  TM  mode  will  see 
a  higher  refractive  index  in  the  poled  regions  than  in  the  un-poled  regions,  becoming 
wavguides  for  TM  light  [61]. 
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The  device  is  now  ready  for  operation.  The  parameters  for  poling  and 
photobleaching  are  different  for  the  various  NLO  polymers;  one  needs  to  determine 
optimum  conditions  through  experimentation.  The  separation  between  the  electrodes  for 
this  example  is  8  (xm.  The  thicknesses  for  the  cladding  and  core  layers  given  above  have 
been  determined,  by  several  researchers,  to  provide  low  propagation  loss  and  are  the 
thicknesses  generally  used  [3, 4, 7,  8, 21, 22, 36, 39-43, 45-47, 48-53, 55-63, 65, 68]. 
However,  an  8  jim  separation  between  electrodes  renders  long  devices  not  suitable  for 
MCM  applications  [35].  This  will  be  addressed  in  a  later. 

5.4.2  Etch  and  Fill 

For  etch  and  fill,  we  once  again  begin  by  depositing  metal  onto  a  substrate  for  the 
ground  plane.  A  5  pm  passive  polymer  cladding  layer  is  spin  cast  over  the  metal  and  heated 
until  dry.  Next,  the  waveguide  circuit  is  patterned  onto  the  cladding  layer  using  standard  IC 
tools  and  techniques.  2  jxm  of  the  cladding  layer  are  etched  using  either  wet  or  dry  etch 
processes.  The  NLO  polymer  is  spin  cast  onto  the  circuit  filling  the  etched  regions  and 
cured  by  heating.  A  3  pm  passive  polymer  cladding  layer  is  spin  cast  on  top  of  the  NLO 
polymer  core  layer  and  cured.  Metal  is  deposited  and  patterned  to  provide  electrodes  over 
the  switching  regions.  The  device  is  then  poled.  The  difference  between  the  two  techniques 
is  that  the  refractive  indices  of  the  passive  polymer  side  cladding  layers  cannot  be  controlled 
as  precisely  as  the  refractive  indices  of  the  photobleached  NLO  polymer  side  cladding 
layers.  This  can  be  important  when  a  precise  interaction  length  L  is  desired.  However,  the 
etch  and  fill  method  does  eliminate  a  metalization,  patterning  and  exposure  step  required  for 
the  photobleach  method. 

5.4.3  Corona  Poling 

A  second  poling  method  commonly  used  is  known  as  corona  poling  [7,  8, 39-41, 45- 
47, 49-54].  After  the  spin  cast  films  are  vacuum  dried  and  heated  to  150°C,  a  small  probe 
tip  emitting  a  high  electric  field,  typically  about  5  KVDC,  is  positioned  2.0  cm  over  the 
heated  NLO  polymer.  This  provides  the  150  VDC/pm  poling  voltage.  The  150°C 
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temperature  and  5  KVDC  poling  voltage  are  maintained  for  2  hours.  The  films  are  cooled 
to  room  temperature,  with  the  electric  field  still  applied.  When  the  films  reach  room 
temperature,  the  voltage  can  be  turned  off,  and  the  material  is  ready  for  EO  device 
fabrication.  The  switching  electrodes  then  need  to  be  deposited  and  patterned.  Corona 
poling  is  also  attractive  for  material  characterization  because  the  NLO  polymer  can  be  poled 
over  a  large  area.  For  device  fabrication,  however,  care  must  be  taken  to  deposit  the  NLO 
material  only  in  the  switching  region,  when  using  corona  poling.  Since  a  poled  NLO 
polymer  exhibits  higher  propagation  loss  than  the  un-poled  NLO  polymer  [3, 4,  8, 21, 22, 
36],  an  EO  circuit  fabricated  using  this  method  will  exhibit  higher  loss  if  the  NLO  material 
is  deposited  anywhere  besides  the  switching  region. 
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5.5  PRACTICAL  NONLINEAR  OPTICAL  POLYMER  ELECTRO-OPTIC 
DEVICES 

Electro-optic  modulators  and  directional  couplers  using  NLO  polymers  have  been 
modelled,  fabricated  and/or  analyzed  extensively  over  the  years  by  many  researchers  [4, 7, 
8, 21, 22, 35, 36,  39-43, 45-70].  The  EO  coefficients  for  various  NLO  polymers  have  been 
measured  by  several  researchers  as  well  [3, 4, 21, 22, 41-43, 45-70, 135]. 

However,  even  with  the  excellent  material  characteristics  that  have  been 
demonstrated,  several  barriers  have  prevented  NLO  polymers  from  progressing  much 
further  than  being  used  for  research  devices.  One  of  these  barriers  is  poling  voltage.  In 
order  to  align  the  moieties  to  attain  the  nonlinearity  of  the  material,  one  needs  to  pole  the 
polymer  with  a  high  voltage.  As  stated  earlier,  the  level  of  this  poling  voltage  can  be  1200 
VDC  or  greater.  Voltage  levels  of  these  magnitudes  prevent  the  ease  of  integration  of  optics 
with  electronics.  Lift-off  techniques  are  being  pursued  to  allow  for  integration,  but  such 
techniques  add  to  the  fabrication  cost  and  can  defeat  the  simplicity  claimed  for  NLO 
polymer  device  fabrication.  An  electronics  company  that  does  not  utilize  lift-off  in  their 
current  fabrication,  might  be  reluctant  to  introduce  it  into  their  process,  unless  there  were  a 
significant  and  marketable  performance  enhancement. 

For  compatibility  with  electronic  processes,  including  solder  baths,  NLO  polymers 
need  to  be  able  to  retain  their  nonlinearity  at  high  temperatures.  It  has  been  demonstrated 
that  the  higher  the  poling  retention  temperature  stability,  the  lower  the  nonlinearity  of  the 
material  [3, 4, 22].  However,  one  would  like  to  achieve  a  high  nonlinearity  for  the  material; 
this  is  its  most  attractive  characteristic.  Due  to  high  poling  voltages  and  lower  nonlinearities 
caused  by  rugedization,  one  should  consider  using  these  materials  in  stand  alone  EO 
devices.  There  is  no  immediate  requirement  to  monolithically  integrate  with  electronics 
and/or  electronic  fabrication  processes.  This  will  allow  for  the  relaxation  of  the  high 
temperature  stability  requirement,  and  it  will  provide  a  polymer  with  a  higher  EO 
coefficient.  It  is  the  belief  of  this  author  that  for  the  near  term,  it  is  more  important  to  get  a 
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competitive  NLO  polymer  device  into  the  market  than  to  focus  on  integration  with 
electronics. 


5.5.1  Electrode  Separation 


Substituting  the  properties  of  the  various  competing  materials,  for  example,  into  Eq. 
(89)  for  the  zero-gap  directional  coupler  (ZGDC)  gives  one  a  good  material  comparison. 
These  are  presented  in  Table  2.  Looking  at  Table  2,  it  is  clear  that  the  NLO  polymer  is  not 
being  utilized  to  its  full  potential.  If  one  could  find  a  way  to  reduce  the  electrode  separation 
for  NLO  polymer  devices,  the  switching  voltage  could  be  reduced.  For  d  =  1  pm,  using  the 
above  parameters,  L  =  1 .5  mm  and  the  voltage-length  is  reduced  to  7.5  V-mm.  Thus,  a 
ZGDC,  using  the  NLO  polymers  that  are  currently  available,  has  the  potential  of  having  the 
lowest  switching  voltage  of  any  integrated  optic  EO  device  demonstrated  to  date.  Table  3 
illustrates  the  potential  of  decreasing  the  electrode  separation.  Looking  at  Table  3,  it  can  be 
seen  that  the  potential  exists  to  realize  extremely  short  devices. 


Table  2.  NLO  Materials  Comparison  [4,  7,  8,  12,  32-35,  37-42,  45, 48,  50,  64,  65, 
136-150], 


Voltage- 

Material  Vn  X  n  r33  d  L  (mm)  Length 


Bulk  m-v 
Semiconductor 

5  V 

830  nm 

3.42 

1.5  pm/V 

2  pm 

5.6  mm 

28  V-mm 

LiNb03 

5  V 

830  nm 

2.2 

32  pm /V 

9  pm 

4.4  mm 

22  V-mm 

MQWm-V 

Semiconductor 

5  V 

830  nm 

3.42 

4.7  pm/V 

2  pm 

1.8  mm 

9  V-mm 

NLO  Polymer 

5  V 

830  nm 

1.63 

26  pm/V 

8  pm 

12  mm 

60  V-mm 
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Table  3.  NLO  Polymer  Comparison  [35] 


Tolerable 

Loss 

V* 

X 

n 

T33 

d 

L  (mm) 

Voltage- 

Length 

12  dB/cm 

5  V 

830  nm 

1.63 

26  pm/V 

1  pm 

1.5  mm 

7.5  V-mm 

21  dB/cm 

5  V 

830  nm 

1.63 

45  pm/V 

1  pm 

0.85  mm 

4.3  V-mm 

187  dB/cm 

5  V 

830  nm 

1.63 

400  pm/V 

1  pm 

96  pm 

0.5  V-mm 

As  a  first  option,  one  could  decrease  the  core  thickness  from  2  pm  to  1  pm.  Doing 
this  would  allow  for  a  1  pm  reduction  in  the  separation  between  electrodes.  However,  this 
would  only  offer  a  slight  reduction  in  the  switching  voltage.  In  order  to  reduce  the  electrode 
separation  any  further,  one  would  have  to  either  reduce  the  thickness  of  the  cladding  layers 
or,  as  in  our  proposal,  render  the  cladding  layers  conductive.  Several  researchers  have 
recently  increased  the  conductivity  of  the  cladding  layers  in  order  to  pole  the  core  layer  in 
the  triple  stack  polymer  device  configuration  to  its  full  NLO  value  [4,  36, 40, 43, 48, 59,  60, 
62],  This  level  of  conductivity  is  not  adequate  enough  to  allow  the  cladding  layers  to  act  as 
switching  electrodes.  By  using  conductive  polymers,  on  the  other  hand,  such  as  HC1  doped 
polyanilene,  AsF5  doped  polyacetylene  or  doped  polypyrrole,  conductivities  of  2xl02  S/cm 
to  lxlO5  S/cm  can  be  realized  [130-134],  The  conductivity  of  polyanilene  is  within  a  factor 
of  300  of  copper.  The  conductivity  of  polyacetylene  is  within  a  factor  of  6  of  copper. 

Increasing  the  conductivity  of  the  cladding  layers  to  these  conductivity  levels, 
however,  increases  the  dielectric  constant  of  the  material.  This  may  produce  significant 
propagation  loss  [119, 120].  It  has  been  demonstrated,  at  X  =  632  nm,  that  undoped,  non- 
conductive  polyanilene,  has  a  dielectric  constant  e  =  8.47  +  1.14i  and  conductive,  HC1 
doped  polyanilene  has  a  dielectric  constant  £  =  13.4  +  2.85i  [120].  The  loss  incurred  by 
rendering  the  cladding  material  highly  conductive  may  prove  too  high  for  practical  device 
operation.  One  might,  therefore,  consider  a  thin  charge  sheet  (<  0.1  pm)  of  conductive 
polymer  spin  casted  between  the  core  and  cladding  layers.  This  would  undoubtedly  increase 
propagation  loss,  but,  not,  one  would  hope,  to  the  levels  anticipated  from  using  conductive 
cladding  layers.  Fig.  20  is  a  schematic  illustrating  the  stacked  structure  under  consideration. 
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Figure  20.  Schematic  of  Triple  Stack  Including  Conductive  Polymer  Charge  Sheets  [35], 

5.5.2  Passive  Zero-Gap  Directional  Coupler 

We  will  first  analyze  our  scheme  for  the  passive  device  (no  voltage  applied).  In 
anticipation  of  the  higher  losses  due  to  the  conductive  polymer,  we  would  only  want  to  use 
the  NLO  polymer  for  the  interaction  region  of  the  switch.  The  input/output  channels,  as 
well  as  the  general  routing,  should  be  fabricated  using  a  passive  polymer,  such  as  General 
Electric  polyetherimide  (ULTEM).  It  has  demonstrated  losses  of  0.24  dB/cm  at  X  =  830  nm 
[2, 23],  Limiting  the  NLO  polymer  to  the  interaction  channel  of  the  device  only  allows  one 
to  relax  the  current  low  loss  requirement  for  NLO  polymers.  Fig.  21  is  a  schematic  of  the 
passive  zero-gap  directional  coupler  under  investigation,  with  the  active  polymer  making  up 
the  interaction  region  and  a  passive  polymer  making  up  the  input/output  channels.  By 
using  the  waveguide  and  material  parameters  given  in  Table  2  and  substituting  d  =  1  pm,  the 
predicted  optimum  interaction  length  (L)  for  the  passive  ZGDC  is  only  101  pm.  Fig.  22  is  a 
BeamPROP  beam  propagation  model  representation  of  this  device. 
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NLO  Polymer 


Passive  Polymer 


Figure  21.  Schematic  of  a  Passive  Zero-Gap  Directional  Coupler  With  an  NLO  Polymer 
Interaction  Region  and  Passive  Polymer  Input/Output  Channels  [35]. 

In  addition,  the  accuracy  of  the  width  of  the  channels  for  a  ZGDC  is  not  too  critical. 
If  we  can  control  the  accuracy  of  the  width  of  the  interaction  channel  to  within  0.1  pm, 
which  is  the  current  tolerance  for  standard  IC  fabrication  tools,  L  will  only  change  by  7.6 
pm.  We  would  be  able  to  control  the  coupling  to  within  93%  of  complete  coupling.  This 
would  prove  very  fabrication  tolerant. 

For  our  ZGDC,  we  are  also  assuming  that  we  will  be  able  to  match  the  refractive 
indices  of  the  NLO  polymer  and  the  passive  polymer  to  within  0.01 .  Calculating  the 
percentage  of  coupling  with  this  amount  of  index  tolerance  predicts  that  we  could  control 
the  coupling  to  within  93%  of  complete  coupling.  This,  again,  proves  very  fabrication 
tolerant.  In  order  to  implement  this  type  of  structure,  we  would  have  to  utilize  the  etch  and 
fill  technique.  However,  with  such  a  fabrication  tolerant  device,  one  does  not  require  the 
precise  refractive  index  control  that  can  be  achieved  using  photobleach.  Fig.  23  is  a  BPM  of 
the  device  behavior,  including  the  0.01  refractive  index  mismatch.  Assuming  1.8  dB  total 
loss  for  this  device,  corresponding  to  1.5  dB/cm  propagation  loss  for  the  conventional 
12  mm  long  device  allows  one  to  tolerate  a  propagation  loss  of  178  dB/cm  for  the  101  jxm 
passive  ZGDC. 
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Transverse  Direction  (£m) 


Figure  22.  Beam  Propagation  Model  of  Passive  Zero-Gap  Directional  Coupler  Using  NLO 
and  Passive  Polymers  With  Matching  Indices. 
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Propagation  Direction  (fim) 
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Figure  23.  Beam  Propagation  Model  of  Passive  Zero-Gap  Directional  Coupler  Using 
NLO  and  Passive  Polymers  With  0.01  Mismatched  Indices. 

5.5.3  Active  Zero-Gap  Directional  Coupler 

For  switching,  we  need  to  implement  the  active  device  (voltage  applied).  In  order  to 
keep  the  switching  voltage  to  a  minimum,  the  active  ZGDC  will  require  a  longer  interaction 
length  than  the  passive  device.  The  coupling  will  also  be  harder  to  predict,  since  any  error 
in  index  matching  or  width  accuracy  will  be  multiplied  by  a  factor  determined  by  dividing 
the  realized  active  device  length  by  the  optimum  passive  device  length.  For  the  active 
device  case,  when  using  the  parameters  from  Table  2  along  with  d  =  1  pm  and  a  0.01 
refractive  index  difference  between  the  NLO  and  passive  material,  our  coupling  error  is  still 
very  low.  Fig.  24  is  a  BPM  of  the  behavior  of  the  active  ZGDC  with  perfect  index  match. 
Fig.  25  is  a  BPM  of  the  active  ZGDC  with  0.01  index  mismatch. 
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Transverse  Direction  (jim) 


Figure  24.  Beam  Propagation  Model  of  Active  Zero-Gap  Directional  Coupler  Using  NLO 
and  Passive  Polymers  With  Matching  Indices. 
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Transverse  Direction  (jim) 


Figure  25.  Beam  Propagation  Model  of  Active  Zero-Gap  Directional  Coupler  Using  NLO 
and  Passive  Polymers  With  0.01  Mismatched  Indices. 

We  would  anticipate  the  coupling  to  be  more  difficult  to  predict  when  fabricating  the 
actual  device,  and  we  would  expect  the  active  device  to  be  less  fabrication  tolerant  than  the 
passive  device.  However,  since  we  will  be  applying  a  voltage  for  switching  anyway,  we 
should  be  able  to  simply  apply  a  voltage  to  null  the  device  in  the  complete  coupling  state.  In 
any  case,  this  1.5  mm  long  device  would  have  the  potential  for  realizing  7.5  V-mm 
switching.  This  is  a  significant  benefit,  which  would  offset  the  inconvenience  of  having  to 
null  the  device.  Low  loss  requirements  for  current  NLO  polymers  could  also  be  relaxed  for 
the  active  ZGDC.  Assuming  a  1.8  dB  total  loss  for  this  device,  the  same  as  we  did  for  the 
passive  ZGDC,  would  allow  us  to  tolerate  a  propagation  loss  of  12  dB/cm  for  the  1.5  mm 
active  ZGDC. 

Hence,  this  length  would  be  short  enough  to  conceivably  build  an  array  NLO 
polymer  EO  devices  within  a  small  chip  for  placement  within  a  multichip  module. 
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5.5.4  Nonlinear  Polymer  to  Passive  Polymer  Coupling 


Another  important  issue  to  be  addressed  is  the  coupling  efficiency  between  the  NLO 
polymer  material  and  the  passive  polymer  material.  If  we  wish  to  use  the  NLO  polymer  for 
the  EO  device  and  the  passive  polymer  for  the  general  routing,  in  order  to  minimize 
propagation  loss,  we  need  to  conceive  better  coupling  techniques  than  the  current  butt 
coupling  technique.  Butt  coupling  two  waveguides  together  has  proven  very  lossy, 
especially  between  two  materials  of  different  refractive  indices.  It  is  not  efficient  enough  for 
coupling  the  NLO  material  with  the  passive  material.  One  possible  method  may  involve 
tapering  the  two  materials  together  or  blending  the  two  materials.  Another  may  be  to  utilize 
the  passive  zero-gap  directional  coupler  or  to  use  grating  couplers.  In  any  event,  in  order  for 
NLO  polymer  devices  to  become  practical  and  competitive  for  computer  applications,  good 
coupling  efficiency  needs  to  be  achieved. 
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6. 


CONCLUSIONS 


Passive  multimode  polymer  waveguides  are  modeled  well  by  the  CAD  program 
BeamPROP.  Active  single-mode  waveguides  do  not  give  as  precise  answers  as  a  Maxwell 
equation  based  three  dimensional  finite  element  method.  However,  when  fabrication 
tolerances  are  taken  into  consideration,  BeamPROP  is  sufficiently  accurate.  Therefore,  in  a 
practical  setting,  BeamPROP  adequately  models  the  behavior  of  the  waveguides. 

Polymer  waveguides  have  a  greater  loss  than  fiber  optical  waveguides.  However,  as 
the  POINT  program  shows,  the  losses  are  low  enough  for  backplane  applications.  POINT 
has  demonstrated  that  multimode  interconnects  are  ready  for  backplane  applications  at  this 
time.  Significant  improvements  are  necessary  for  active  single  mode  waveguides  based  on 
polymers  to  overcome  the  lead  that  lithium  niobate  based  electro-optical  modulator 
waveguides  have  at  this  time. 


A  scheme  for  realizing  1-2  mm  long  EO  devices,  utilizing  currently  available  non¬ 
linear  optic  polymers,  was  introduced.  It  offers  the  potential  to  achieve  the  lowest  switching 
voltage  of  any  integrated  optic  EO  device  of  that  length  demonstrated  to  date.  The  major 
factors  are  the  propagation  loss  dependance,  using  the  conductive  polymer  cladding,  and  the 
coupling  efficiency  between  the  NLO  polymer  and  passive  polymer.  If  losses  are  not 
excessive,  these  devices  have  excellent  potential  for  success.  If  future  EO  coefficients  can 
then  be  pushed  higher,  which  would  relax  the  current  low  loss  constraint  and  the  current 
high  temperature  stability  constraint  even  further,  there  exists  the  potential  for  100  pm  long 
active  EO  devices  operating  at  TTL  voltage  levels  (see  Table  3).  This  would  also  be  a 
fabrication  tolerant  device  that  would  not  require  a  nulling  voltage.  There  are  undoubtedly 
other  possible  methods  for  reducing  the  electrode  separation,  but  the  ultimate  goal  should  be 
to  reduce  that  separation.  This  is  the  only  way  to  produce  practical  and  competitive  electro¬ 
optic  devices  using  available  non-linear  optic  polymers. 
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